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SUMMARY 


FAR-WING  COLLISION  BROADENING  STUDIES  -  Chapter  I 

The  work  reported  in  Chapter  I  is  a  continuation  of  our  studies  of 
absorption  in  the  far  wings  of  molecular  bands.  It  is  in  these  spectral 
regions  that  conventional  pressure  broadening  theories  of  photo  absorption 
are  Known  to  be  inapplicable.  A  striking  example  is  to  be  found  in  the  far 
wings  of  certain  infrared  bands  of  C02.  Here  it  is  found  that  deviations 
from  conventional  (Loreotz-type)  line  shapes  can  assume  order  of  magnitude 
proportions  at  spectral  distances  of  the  order  of  100  to  200  cm"1  from  the 
band  center. 

The  present  work  is  focused  on  the  markedly  sub-Lorentzian  absorption 
which  has  been  observed  in  the  far  wings  of  the  C02  band  centered  at  2350  cm-1 
(4.3  microns).  The  basic  broadening  mechanism  is  assumed  to  arise  from  per¬ 
turbations  of  the  rotational  motion  via  short-range  repulsive  interaction 
between  the  absorbing  molecule  and  a  perturbing  system  (assumed  monatomic  for 
simplicity).  In  the  limit  in  which  the  repulsive  interaction  is  taken  to 
have  infinite  slope  (hard-sphere  limit),  the  absorption  is  found  to  vary  with 
frequency  as  (u>-oj0)"^.  In  the  more  realistic  case  in  which  the  repulsive 
interaction  is  assumed  to  be  of  the  form,  V  =  e~ar  (with  realistic  values  for 
the  "range-constant,"  a),  it  is  found  that,  in  the  far-wing,  the  absorption 
may  be  written  as  the  product  of  a  Lorentz  line-shape  and  a  factor  which  falls 
off  exponentially  with  increasing  frequency-shift  from  the  line  center,  in 
general  agreement  with  Benedict's  empirical  formula.  In  fact,  an  upper- 
bound  estimate  in  the  frequency  range,  w-w0  =  200-250  cm"1,  is  lower  than  the 
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SUMMARY  (continued) 


measured  absorption  by  a  factor  which  varies  from  10  to  15.  It  is  strongly 
believed  that  this  discrepancy  is  due  principally  to  the  omission  of  at¬ 
tractive  interactions;  their  inclusion  into  the  mathematical  apparatus  of 
the  theory  is  feasible,  and  will  constitute  the  subject  matter  of  a  future 
investigation. 

BEAM  DYNAMICS  -  Chapter  II 

We  have  obtained  a  fairly  simple  closed  form  expression  for  the  maximum 
flux  which  the  thermal  self-defocusing  effect  permits  us  to  transmit  through 
the  atmosphere.  This  analysis  takes  into  account  the  effect  of  target  motion, 
and  is  an  extension  of  a  previous  analysis1  in  which  only  stationary  targets 
are  considered.  In  the  case  of  a  stationary  target  the  limitation  was  found 
to  be  on  the  energy  per  unit  area  which  can  be  transmitted,  whereas  for  a 
moving  target  the  limitation  is  found  to  be  on  the  power  density.  The  ex¬ 
pression  for  the  maximum  flux  depends  upon  such  parameters  as  the  initial 
beam  intensity,  the  initial  beam  diameter,  the  target  distance,  the  rotational 
rate  of  the  beam  *nd  the  effective  absorption  coefficient  for  heating  of  air 
by  light  at  the  laser  frequency. 

Our  analysis  is  in  the  ray  optics  approximation,  so  that  diffraction 
spreading  is  Ignored.  Furthermore,  to  avoid  the  complication  of  accounting 
for  sonic  effects,  we  have  assumed  that  the  time  required  for  the  beam  to 
pass  through  any  air  mass  is  longer  than  the  time  required  for  sound  to 
propagate  in  air  across  the  beam;  as  a  consequence,  our  results  are  meaningful 
only  for  targets  moving  at  subsonic  velocities.  In  addition,  rather  than 
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SUMMARY  (continued) 


treat  the  rotation  of  the  beam  exactly  by  using  time-varying  Initial  con¬ 
ditions  on  the  ray  trajectories,  we  have  simulated  It  by  assigning  to  each 
air  mass  an  effective  illumination  period,  which  depends  on  its  distance 
from  the  laser. 

In  general  the  manner  in  which  the  beam  self-defocuses  depends  upon 
the  specific  shape  of  its  initial  power  profile.  We  have  concerned  our¬ 
selves  with  the  special  case  in  which  the  initial  flux  distribution  is 
uniform  in  the  interior  of  the  beam,  and  falls  off  exponentially  at  the  edge, 
as  illustrated  in  Fig.  1.  For  this  case  we  find  that  the  interior  rays 
propagate  undeflected  up  to  a  certain  distance,  then  deflect  outwards;  the 
closer  the  ray  is  to  the  edge,  the  shorter  the  distance  is  at  which  It  first 
deflects  (see  Fig.  2).  Therefore,  the  behavior  of  the  beam  may  be  described 
as  follows:  The  flux  in  the  interior  remains  undiminished,  but  the  outer 
rays  gradually  peel  off,  so  that  the  diameter  of  the  undegraded  interior 
decreases  with  distance. 

Our  solution  is  valid  only  in  certain  regions  of  space;  as  a  result  we 
cannot  describe  the  behavior  of  the  rays  after  they  begin  to  deflect,  but 
we  can  predict  the  diameter  of  the  undegraded  beam  interior  as  a  function  of 
propagation  distance.  Therefore,  as  a  criterion  for  when  the  beam  has  lost 
its  effectiveness,  we  have  assumed  somewhat  arbitrarily  that  the  beam  edge 
may  be  permitted  to  eat  its  way  into  the  Interior  by  no  more  than  one  fourth 
of  the  initial  beam  diameter. 

Using  this  criterion  we  have  provided  several  illustrative  numerical 
examples;  In  all  of  these  we  have  assumed  a  C02  laser  with  a  diameter  of 
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one  meter.  For  a  stationary  target  in  the  zenith  direction  at  an  altitude 

much  greater  than  the  atmospheric  scale  height  (8  km)  we  find  a  limiting 

2 

energy/area  of  0.95  joules/cm  .  For  a  stationary  target  on  the  ground  at 
a  distance  of  1  km  from  the  laser  we  find  a  limit  of  240  joules/cm^.  If 
the  target  is  flying  directly  overhead  at  an  altitude  of  10  km  and  with  a 

O 

velocity  of  500  mph,  we  obtain  a  limiting  flux  of  197  watts/cm  .  Finally, 

if  the  target  is  on  the  ground  and  moving  normal  to  the  line  of  sight  with 

a  velocity  of  10  m/sec  (=  21.6  mph)  at  a  distance  of  1  km,  we  estimate  a 

2 

maximum  flux  of  610  watts/cm  . 
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Chapter  I 

FAR-WING  COLLISION  BROADENING  STUDIES 
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1. 


INTRODUCTION 


In  recent  years,  there  has  developed  an  awareness  of  the  fact  that 
published  theories  of  pressure  broadening  of  spectral  lines  are  quite  at 
variance  with  experiment  with  respect  to  magnitude  and  frequency-dependence 
of  absorption  In  the  extreme  wings  of  the  lines.  In  particular,  In  the  case 
of  Infrared  vibration-rotation  molecular  bands,  the  theories  universally  pre¬ 
dict  the  Lorentz-line  shape;  their  principal  focus  has  to  do  with  the  dependence 
of  the  two  parameters  -  width  and  shift  -  of  the  individual  rotational  com¬ 
ponents  on  such  variables  as  J-values,  temperature,  and  types  of  assumed  Inter- 
molecular  Interactions  Experimentally,  however.  It  is  found  that  deviations 
from  the  Lorentz  line-shape  occur  at  spectral  distances  from  the  center  of  an 
Individual  line  as  small  as  a  few  cm"1.  In  the  extreme  wings  -  at  distances 
100  -200  cm"1  from  the  band  center  of  vibration-rotation  lines  -  the 

discrepancies  may  assume  order-of-magnltude  proportions.  In  fact,  recent 

1  2 

experiments  by  Winters,  Silverman,  and  Benedict  and  Burch  and  Gryvnak  have 
shown  that  the  absorption  in  the  extreme  high-frequency  wings  of  a  number  of 
CO2  bands  is  sub-Lorentzian  by  several  orders  of  magnitude  (i^lO  ).  Benedict 
et  al1  have  in  particular  found  that.  In  order  to  fit  their  experimental  data. 

It  was  necessary  to  multiply  the  Lorentz  shape  by  an  empirical  factor  which 

decreases  exponentially  with  increasing  distance  from  the  band-center.  Similar 

? 

results  have  been  obtained  by  Burch  and  Gryvnak  ,  with,  however,  one  noteworthy 
exception  Namely,  in  the  far-wing  (high-frequency)  wing  of  the  S^-band  of 
COg  In  the  vicinity  of  7000  cm"1,  the  absorption  of  the  self-braodened  band, 
while  sub-Lorentzian,  is  much  higher  than  that  found  by  Benedict  et  al1  for 
the  comparable  situation  of  the  self-broadened  wing  of  the  fundamental  ^  band 
(in  Me  vicinity  of  2400  cm'1). 


During  the  past  year,  work  has  been  focused  primarily  on  developing 
a  basic  understanding  of  the  above  described  sub-Lorentzian  behavior  in  the 
far  wings  of  vibration-rotation  bands.  This  work  constitutes  the  subject 
matter  of  the  present  report 

Two  basic  mechanisms  have  been  considered:  (1)  collisional  pertur¬ 
bation  of  the  rotational  motion  of  the  absorbing  molecule,  and  (2)  collisional 
perturbation  of  its  vibrational  motion  In  our  work,  the  main  effort  has 

3 

been  concentrated  on  the  first  of  these,  namely  rotational  perturbations. 

The  basic  treatment  is  given  in  Section  2.  Comparison  with  experiment  (and 
discussion  thereof)  is  presented  in  Section  3;  in  Section  4  future  possible 
developments  of  the  theory  are  outlined  and  overall  prospects  assessed 

With  respect  to  the  mechanism  of  vibrational  perturbations,  brief 

preliminary  estimates,  based  on  a  simplification  of  the  treatment  given  in 

a 

one  of  our  previous  ONR  reports  indicate  that  its  contribution  to  the  far- 
wmg  absorption  is  substantially  smaller  than  that  due  to  rotational  pertur¬ 
bations  A  systematic  investigat  on  of  this  mechanism  is  therefore  reserved 
for  the  future 
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2.  ABSORPTION  DUE  TO  ROTATIONAL  PERTURBATIONS  -  BASIC  THEORY 

In  this  section,  the  far-wing  absorption  due  to  collisional  perturbations 
of  the  rotational  motion  of  a  radiating  molecule  will  be  studied.  In  line 
with  the  remarks  of  the  introduction,  the  radiating  molecule  will  be  taken  to 
be  a  linear  (COg-lIke)  system;  for  simplicity  the  perturber  will  be  assumed 
monatomic  (noble-gas  atom). 

The  starting  point  of  the  calculation  is  a  vector  generalization  of 

c 

equations  (1)  and  (2)  of  a  previous  paper"  for  the  spectral  intensity,  I(u>); 
it  takes  the  form 


M«)  = 

<JA(«)| 

/Av  • 

(2.1) 

A(u)  * 

1 

a+T/2  . i  \ . 

-i.(w-u)  )t 

e  r^(t)dt  . 

(2.2) 

(2irT)is 

J-T/2 

As  in  Ref.  5  the  symbol  denotes  an  average  over  all  types  of 

collisions,  and  T  is  an  arbitrarily  large  time  interval.  The  new  ingredient 
is  the  replacement  of  the  scalar  phase-shift  factor,  e^^,  by  the  unit 
vector,  njt),  giving  the  direction  of  the  transition  dipole-moment,  ip(t),  of 
the  radiating  molecule  via  the  expression, 

iu  t 

u_(t)  =  n(t)  e  0  ,  (2.3) 

where  is  the  absolute  magnitude  of  the  transition  di pole-moment^  and  where 
uQ  is  the  frequency  at  the  band  origin,  i.e.,  the  frequency  associated  with 
purely  vibrational  motion.  The  time-variation  of  nU)  is  thus  associated 
solely  with  the  rotatory  motion  of  the  axis  of  symmetry  of  the  linear  molecule; 
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in  line  with  the  experiments  discussed  in  the  introduction,  the  treatment  will 
be  confined  to  transitions  involving  the  so-called  v3«mode  (1 . e. ,  asymmetric 
vibrations  parallel  to  the  axis  of  the  linear  molecule),  so  that  n(t)  Is 
parallel  to  the  axis  of  symmetry. 

It  may  be  remarked  that,  with  |n(t)|  *  1,  the  factor  l/^uT)*5  In  (2.2) 
guarantees  that  I(w)  is  properly  normalized,  i.e., 


f°° 

I(u)do)  =  1  .  (2.4) 

o 

in  the  present  treatment,  the  rotational  motion  will  be  treated  clas¬ 
sically;^  one  then  has 


dn(t) 

dt 


u)p(t)  x  n  , 


(2.5) 


where  u)p(t)  is  the  instantaneous  rotational  frequency. 

Performing  two  successive  integrations-by-parts,  and  introducing  (2.5) 
one  has 


A(w)  =  ( 2-n-T ) 


-k 


f+T/2  w_xn(t)  -i(w-u)0)t 

)  e 

T/2  0 


dt 


*  WVh 


+T/2 


J. 


duu(t)  .1 

*  n(t)+aip(t)  x  ij)p(t)xn(t)]J  -i(a)-w0)t 

e 


dt  , 


T/2 


(2.6) 


where  the  integrated  terms  at  T/2  and  -T/2  have  been  discarded. 
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The  basic  motivation  in  going  from  (2.2)  to  (2.6)  is  contained  in  the 
classical  picture  of  rotational  perturbations,  a  picture  originally  developed 
by  Debye,  and  which  we  feel  particularly  appropriate  for  the  description  of 

Q 

perturbations  arising  from  short-range  (repulsive)  interaction  potentials 
between  the  radiating  molecule  and  a  (monatomic)  perturber.  According  to 
this  picture,  the  effect  of  such  collisions  is  simply  to  change  the  vector 
angular  velocity  of  rotation,  (generally,  both  in  direction  and  magnitude). 
Following  Debye,  we  shall  assume  (at  least  for  the  time  being)  that  such 
changes  occur  instantaneously  (impact-type  collisions);  subsequently,  a 
generalization  taking  into  account  the  finite  time-duration  of  collisions 
will  be  introduced. 

Accordingly,  let  it  be  assumed  that  the  molecule  undergoes  a  series  of 

collisions  at  times  t^ (-T/2  <  tj  <  t^ <  .  <  t^  <  .....  <  tR  <  T/2),  at 

which  a  changes  suddenly.  In  particular,  the  angular  velocities  before  and 

r 

after  the  i'th  collision  are  respectively  denoted  as  -  Au^/2  and 

+  %j/  2;  thus,  Awj  is  the  change  at  the  i'th  collision  and  is  the 
arithmetic  mean  of  the  angular  velocities  before  and  after  the  i'th  collision. 
It  is  then  immediately  apparent  that,  for  the  factor  da^/dt,  occurring  in  the 
first  tenri  of  the  curly  bracket  of  (2.6),  one  may  write 

du  N 

~  1  Aw.j  6  ( t-t  ■■ )  ,  (2.7) 

with  the  result  that  the  contribution,  A^U),  of  this  first  term  to  A(w)  is 


A 


(I) 


U) 


N 

■  l  (2rT)-* 
i  =  l 


Aolk  xn^j 

U-^0)2 


-i  (u>-u>0)t^ 


(2.8) 
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where  Is  the  (momentary)  orientation  of  the  axis  of  the  linear  molecule  at 
the  time  of  the  I'th  impact,  t^. 

In  order  to  treat  the  second  member  of  (2.6)  in  this  spirit,  let  us 
carry  out  another  integration-by-parts  on  it.  Discarding,  again,  the  physically 
spurious  contributions  at  ±  T/2,  we  have 

tOpXn(t)]  +  ojx 
i(ui-u>0) 

+  Rem  (2.9) 

where  the  term  denoted  as  "Rem"  is  obtained  by  differentiating  ni(t)  and  using 

(2.5);  its  integrand  is  a  vector-product  expression  involving  the  factor, 

(oj— u)Q )  ,  three  powers  of  w_(t),  the  orientation-vector  njt),  and  the  ex- 
-i(u)-u)  )t 

ponential  e 

It  is  clear  from  these  remarks  that  "Rem"  may  be  treated  by  an 
integration-by-parts  to  yield  terms  whose  ’ntegrands  have  magnitudes 
together  with  "remainder-integrands"  ^  oj£/(u)-u>0)  .  Iterating  this 
procedure,  one  may  in  fact  construct  a  series  in  ascending  powers  of  u>p(t)/(u)-u;0) 
it  is  here  assumed  without  proof  that  such  a  series  is  at  least  asymptotically 
convergent  (i.e.,  in  that  the  error  is  of  the  order  of  the  first  discarded  term). 
The  applicability  of  the  whole  procedure  then  depends  upon  the  smallness  of  the 
parameter,  u)p/(u-w0),  where  uip  is  a  typical  rotational  frequency.  Assuming 
that  this  parameter  is  sufficiently  small  let  us  discard  "Rem"  in  (2.9)  without 
further  ado. 


additional 
^  MPdtuP/dt 
(w-w0)4 
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The  first  term  on  the  r.h.s.  of  (2.9),  herewith  designated  as  A^(w), 
Is  now  evaluated  straightforwardly  with  the  use  of  (2.7);  one  has 


A(2)(u,) 


(2.10) 


Equation  (2.10)  Is  to  be  regarded  as  a  first-order  correction  to  A^(w)  in 
the  (assumed)  asymptotic  series 

AU)  «  A(1)U)  +  A(2)(u>)  + .  (2,11) 


which  has  now  to  be  inserted  into  (2.1). 

In  carrying  out  this  last  mentioned  stage  of  the  calculation,  one  may 
achieve  a  substantial  simplification  by  introducing  the  so-called  "wing- 
assumption,"  specifically,  that 


(ui-U>0)  T  »  1 


(2.12) 


where  t  is  the  mean  time  between  successive  collisions  (assumed  to  be  binary) 

-i  (u>-u)  )t.j 

This  condition  guarantees  that  the  phase  factors,  e  (after  averaging 

over  the  appropriate  statistical  distribution  of  times  between  collisions, 

*1+1  "^i*  are  essentia11y  uncorrelated,  i.e.. 


(2.U) 


Ignoring  these  correlations,  It  then  follows  (subject  to  one  further  qualifi¬ 
cation,  to  be  discussed  immediately  below)  that  cross-terms  in  (2. 1 )(invol ving 
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products  of  contributions  of  different  collisions)  drop  out,  and  one  has 


I(«) 


Ao>j  x  n4 

— — -7  +  ( 

(oi-U)0) 


2 

Av 


(2.14) 


where  ( . )i  represents  the  i'th  term  In  the  r.h.s.  of  (2.10)*®. 

The  above  mentioned  qualification  concerns  the  possibility  that  the 
orientation  vector,  njt),  may  undergo  variations  between  collisions  which 

i  ( ijO—tO  )  t 

compensate  those  of  the  phase  factor,  e  0  .  It  is  however,  immediately 
seen  from  (2.5)  that  this  possibility  can  only  be  realized  provided  that 
there  is  available  a  rotational  frequency,  u)p,  in  the  near  vicinity  of  (u>-u>0), 
such  that 

(u>-u0  -  U>p)t  1  . 

But  in  this  case  our  spectral  frequency  would  be  located  close  to  one  of  the 
unperturbed  rotational  components  of  the  band  -  a  situation  which  is  anti¬ 
thetical  to  that  of  interest  in  the  present  paper. 

Returning  to  (2.14),  one  notes  that,  since  the  collisions  are  statis¬ 
tically  independent  of  each  other,  the  sum  may  simply  be  replaced  by  the 
total  number  of  collisions  in  time  T;  since  this  number  N  is  T/t,  one  has 


K«) 


Au)j  x  n4 

-4*  ( 

(u-u>0) 


f 


(2.15) 
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where  the  term  ( . has  been  discarded  In  the  Interest  of  obtaining  the 

simplest  zero'th  order  result.**  With  perpendicular  to  (as  In  the  case 
for  all  angular  velocities  of  rotation  of  linear  molecules),  (2.15)  may  be 
written  as 


At  this  point,  some  brief  remarks  on  the  quantity,  1/t,  are  In  order. 

As  remarked  above  (and  as  will  be  developed  much  more  fully  below),  an  impact- 
type  approach  is  only  valid  when  the  duration,  tc,  of  a  collision  Is  small 
compared  to  a  time  which  (cf.  especially  (2.35)  and  (2.36)  and  the  text  para¬ 
graph  after  that  containing  these  equations)  Is  ^  — — r  ,  i.e., 

|(D-U>0| 

Iu-WqIt^  <<  1  .  (2.17) 

Now  It  turns  out  (as  will  be  seen  below)  that  for  values  of  |w-ui0|  of  ex¬ 
perimental  relevance  (^  200  cm"*),  (2.17)  can  be  satisfied  only  for  collisions 
governed  by  steeply  rising,  short-range  repulsive  potentials  (in  fact,  rather 
more  steeply  rising  than  is  realistically  to  be  expected).  Moreover,  It  will 
also  be  found  that,  when  (2.17)  Is  not  satisfied,  the  predictions  of  Impact 
theory,  as  given  by  (2.16),  are  modified  principally  by  the  Inclusion  of  an 
exponential  factor  whose  argument  Is  ^  2*|u-w0|tc. 

From  these  remarks.  It  follows  a  fortiori  that  collisions  of  Impact 
parameter  appreciably  larger  than  gas-kinetic  radii  -  I.e.,  those  Involving 
the  multi  pole- type,  long  range  Interactions  (of  predominant  Interest  to  the 
Lorentz-type  Impact  theories,  whose  domain  of  validity  Is  confined 
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to  the  near  vicinity  of  an  Individual  rotational  component  of  a  band  (<  1  cm”*) 
-will  be  quite  Ineffective  In  the  far-wlng  region.  Thus,  in  expressing  1/t 
In  terms  of  a  cross-section,  a  ,  via  the  relationship 

i  =  Nvo  ,  (2.18) 

(where  v  Is  a  suitably  averaged  velocity  of  the  perturbing  atom  relative  to 
the  radiating  molecule),  we  must  Infer  from  the  above  remarks  that  a  Is  at 
most  of  the  order  of  a  gas-kinetic  cross-section. 

Excursus:  At  this  point.  It  Is  of  Interest,  and  perhaps  useful,  to  discuss 

12 

certain  other  features  of  comparison  between  the  Lorentz-type  Impact  theories 

and  that  contained  In  Eq.  (2.16).  As  stated  above  (cf.  Section  1),  and  ex- 

12 

pllcltly  In  the  cited  references,  the  Lorentz-type  theories  all  predict 
line  shapes  whose  wing-behavior  Is  of  the  form 


I(u)  « 


(2.19) 


Now  a  fundamental  postulate  of  these  theories  Is  that  collisions  be  regarded 
as  Impacts  -  l.e.,  of  time  duration  small  compared  to  any  other  physically 
relevant  time- like  quantity  In  the  problem.  Our  theory  -  or  at  least  the 
version  presented  up  to  this  point  -  Is  based  on  the  same  postulate.  Never¬ 
theless,  as  Is  clearly  seen  from  (2.16),  the  frequency-dependence  Is  of  the 


type 


Wherein  lies  the  difference? 


The  answer,  In  our  opinion.  Is  contained  In  an  additional  assumption, 

13 

first  Introduced  Implicitly  by  Anderson;  namely,  if  one  Is  dealing  with  a 
number  of  lines  (or  individual  components  of  a  composite  line)  of  frequencies 
wj...  It  Is  necessary  that  the  spectral  frequency,  w,  lie  closer  to  any 
one  of  the  lines,  say  <^,  than  to  the  others,  l.e.. 


<<  •  (2.21) 
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In  our  theory,  just  the  opposite  assumption  is  made,  namely 

Ju-u»0l  >>  l&u^l 


It  Is  hoped  to  verify  the  above  conjecture  In  detail  in  the  future. 

At  this  time  we  proceed  to  discuss  matters  rather  more  relevant  for  the  Inter¬ 
pretation  of  the  C02  wing-experiments.  Primary  among  these  Is  the  question 
of  the  validity  of  an  Impact-type  theory  (as  described  by  Eq.  (2.7)  and  the 
two  paragraphs  preceding  It).  Specifically,  we  now  consider  the  effect  of  the 
actual  time-duration  of  collisions  -  l.e.,  the  time  Interval  over  which  the 
collision- Induced  change  In  rotational  frequency,  Ak>jr  occurs. 

In  order  to  Investigate  this  question  it  Is  necessary  to  determine  the 
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details  of  the  time  variation  of  Up(t)  in  the  Interval  of  collision.  To  do 
this  properly  would  involve  rather  complicated  calculations  and  accuiate 
knowledge  of  Interatomic  interactions  (between  the  perturbing  atom  and  trose 
of  the  radiating  linear  molecule).  Since  the  latter  Information  Is  not  avail¬ 
able,  and  since  we  wish  to  get  at  the  basic  essence  of  the  effect,  we  shall 
attack  the  problem  In  an  admittedly  simple-minded  manner,  as  follows. 
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Let  us  begin  by  noting  that,  as  remarked  in  the  text  surrounding 
Eq.  (2.17),  we  have  already  restricted  the  treatment  to  a  model  In  which  the 
Interaction  V1nt,  between  the  radiating  molecule  and  perturbing  (nohle-gas) 
atom  Is  repulsive.  In  the  Interest  of  definiteness  and  simplicity,  we  now 
postulate  a  particular  form  for  this  Interaction.  SpecIflcaHy,  It  Is  as¬ 
sumed  that  V.  *  Is  the  sum  of  Interactions  between  the  noble-gas  perturber 
and  the  Individual  atoms  of  the  radiating  molecule,  i.e., 

»,*  *  \  »u  ''S>r,a  •  (2-23) 

(where  r.  Is  the  Interatomic  distance  between  the  perturber  (sympollzeo  by 
subscript  "a")  and  the  I'th  atomic  constituent  of  the  molecule,  and  the  Via 
and  a.  are  constants  determining  the  magnitudes  and  ranges  of  the  component 

I  a 

Interactions!. 

We  have  now  to  note  a  point  of  decisive  Importance  for  the  calculation. 
Namely,  the  range  constants,  a^a,  turn  out  to  be  generally  of  the  order  of 
4xl0®cm“ *  (I.e.,  a  few  Inverse  angstroms),  whereas  a  typical  interaction 
distance,  Ria,  defined,  e.g.  in  terms  of  the  distance  of  closest  approach  of 
two  colliding  atoms)  Is  v  3xl0"8  i.e.,  a  few  angstroms;  one  thus  has 


°1aR1a 


1 

T7 


(2.24) 


The  smallness  of  the  parameter,  lA»iaRia.  permits  certain  substantial 
simplifications.  The  most  Important  of  these  is  associated  with  the  fact 
that  she  depth  of  penetration  of  two  colliding  atoms,  which  Is  *-  l/»1a. 
which,  as  will  be  seen  shortly.  Is  intimately  linked  up  with  the  detailed 
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time-dependence  of  In  a  collision  Interval,  Is  small  compared  to  the  radius 
of  curvature  of  surfaces  of  constant  potential.  This  feature  suggests  that  a 
reasonable  approximation  would  be  to  neglect  said  curvature  altogether.  In 
particular.  It  suggests  that  the  time-dependence  of  the  relative  velocity  of 
the  noble-gas  atom  with  respect  to  Its  collision  partner1®  (e.g.,  In  the  case 
of  COg,  one  of  the  oxygen  atoms)  be  treated  as  a  one-dimensional  problem  of 
a  particle  moving  In  a  potential  field  of  the  form  V  e"aX  (where  the  direction 
of  the  one-dimensional  coordinate,  x.  Is  taken  parallel  to  the  line  joining 
the  centers  of  the  two  colliding  atoms).  For  such  a  problem,  It  Is  Inmedlately 
seen  that  only  the  normal  component  of  relative  velocity  (parallel  to  the 
x-axIs)  would  be  altered  by  collision;  denoting  this  component  simply  by  v, 
one  has 


where  xc  Is  the  distance  of  closest  approach,17  and  vQ  the  normal  component 
of  relative  velocity  at  Infinity  (l.e.,  outside  the  penetration  region).  The 
solution  of  (2.25)  Is 


which  gives 


2  aVo^t'tc^ 

*  xc  +  q  lo9  cosh  — T - 


v  •  vQ  tanh  — g -  , 


(2.26) 


(tc  being  the  time  at  which  the  distance  of  closest  approach  Is  attained). 
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Before  proceeding  further,  we  must  point  out  that  the  relevance  of  the 
simple  one-dimensional  problem,  discussed  Inmedlately  above,  to  our  problem 
rests  upon  an  additional  approximation  In  which  the  repulsive  potential  field 
Is  regarded  as  stationary  (in  the  relative  coordinate  of  the  two  atoms)  during 
the  time  Interval  of  a  collision  ('v  l/avQ).  Such  an  Idealization  would  be 
valid  provided  that,  apart  from  their  Interaction,  the  two  colliding  atoms 
could  be  considered  as  free,  and  hence  moving  In  straight  lines.  Actually, 
however,  one  of  the  atoms  Is  also  rotating;  taking  this  complication  Into 
account  would  require  us  to  consider  a  problem  In  which  the  (planar)  potential 
energy  surface  is  undergoing  rotatory  motion.  The  neglect  of  such  motion 
requires  that  the  parameter  uip/ov Q  be  sufficiently  small,  l.e.. 


(i) 

P 

aVo 


<< 


1  . 


(2.27) 


It  has  now  to  be  remarked  that  (2.27),  which  Is  equivalent  to  requiring 
that  the  angle  of  rotation  during  a  collision  interval  be  small  (compared  to 
a  radian),  may  be  somewhat  more  marginally  satisfied  In  practice  than  (2.24). 
For  example.  In  Cl^,  «  vrot^RC0  where  the  ^  distance  h  1  angstrom) 

In  this  case,  (2.27)  becomes 


vrot  1 


1 

T 


vrot/vo 


1  . 


How,  In  the  worst  possible  case,  in  which  the  noble-gas  atom  is  very  heavy, 
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vrot  *  v  and  one  has  v^/  ^  1  /4,  wMch  is  less  satisfactory  than  (2.24). 

However,  in  a  treatment  of  the  type  given  In  this  paper,  it  appears  to  us  as 
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still  appropriate  to  regard  u  /av  as  a  small  parameter;  higher  order  corrections 

r 
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can  always  be  studied  by  perturbation  techniques  at  some  future  occasion. 

Within  the  framework  of  conditions  (2.24)  and  (2.27),  we  now  postulate 
that  the  time-dependence  of  u)p( t)  in  the  interval  of  a  collision  is  the  same 
as  that  given  by  the  r.h.s.  of  (2.26),  i.e., 

Au.  av . (t-t. ) 

Wp{t)  =  Wpi  +  tanh  — ,  (2.28) 


where  v,  is  the  value  of  v  for  the  i'th  collision. 

1  o 

The  basis  of  this  postulate  is  contained  in  the  following  physical 
argument.  Let  us  consider,  for  example,  the  collision  of  a  noble-gas  atom 
with  one  of  the  oxygen  atoms  of  a  C02  molecule.  Equation  (2.26)  implies  that, 
during  the  collision-interval  there  exists  a  force,  _F,  given  by 


F 


fo  sech 


«Vi(t-ti) 


(2.29) 


where  [within  the  Tramework  of  (2.24)  and  (2.27)],  is  a  temporally-constant 
vector  parallel  to  the  line  of  centers  of  the  two  colliding  atoms  (at  the 
time  instant,  ti ,  of  closest  approach  in  the  i'th  collision).  Now,  the 
component  of  this  force  perpendicular  to  the  C02  axis  gives  rise  to  a  torque, 
which  [again  within  the  framework  of  (2.24)  and  (2.27)]  is  of  the  form 


T 


T  sech 
-o 


aV.(t-ti) 
- 2 - 


(2.30a) 


wnere 


lo 


=  Rco  £  x  £o 


(2.30b) 
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(2.31) 


(with  nj  the  value  of  £  at  t  =  t^).  Noting,  now,  that 


=  T 


(with  I  =  moment  of  inertia  of  the  molecule),  one  integrates,  obtaining 
(2.28)20,  q.e.d. 

The  stage  has  been  reached  where  our  program  of  going  beyond  the  confines 

of  impact  theory  may  now  be  realized.  Let  us  first  obtain  a  lowest-order 

result  by  substituting  (2.28)  into  the  first  term  of  (2.6).  One  has  then 
21 

to  evaluate 


-Aj(u)  = 


(2*T )-* 


+T/2 


-T/2 


o V.(Au).X  n.)  ,  ov.  (t-t.)  -i(a)-!D„)t 

--1  ~4  _1  seclv  e  0  dt. 


(2.32) 


2 

Noting  that  the  sech  factor  in  (2.32)  limits  the  time  integration  to 

intervals,  ^  2/av.,  which  are  small  compared  to  the  time  between  collisions, 

i  (w-oj  )t 

and  that  the  phase  factor  e  undergoes  many  oscillations  between 

collisions,  one  has 


-Ai(u>) 


Au^  x  n_.. 


where 


aVi 

T* 


.+00 


sech 


ov^t 


-i (w-w  )t 

e  0  dt  , 


(2.33) 
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Is  evaluated  in  Appendix  I;  the  result  [cf.  (1.5)1 


o  /  .  -ir(w-u)  )/aV4 

2*(a)-w0)  fi  0"  1 

oTj  -2ir(ui-aJo)/av^ 


(2.34) 


The  contribution  of  AjU)  to  I(u)  is  obtained  by  the  same  procedure  as 
was  used  earlier  [cf.  text  between  Eqs.  (2.11)  and  (2.15)];  the  result  is 


(2.35)  embodies  two  approximations.  The  first  of  these  is  the  neglect  of 

the  term  ( . )^  [which  represents  the  contribution  of  the  second  term  of 

(2.6)1.  Unlike  the  impact  case,  discussed  above,  the  smallness  parameter  is 
not  /u),  but  rather  law  |av4;  the  second  is  based  in  the  smallness  of  the 

P  P  I 

parameter  ^../av^,  in  thet  in  the  first  term  of  (2.6)  n  is  taken  to  be 
constant  (equal  to  its  value,  n^  at  the  instant  of  closest  approach  of  the 
colliding  atoms  (cf.  text  surrounding  (2.27)  and  following  (2.28),  especially 
footnote  19. 
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Equation  (2.35),  in  conjunction  with  (2.36),  contains  the  looked-for 
generalization  of  (2.15),  i.e.,  the  inclusion  of  the  effect  of  the  finite 
duration  of  collisions,  tc  ^  2/av^  (cf.  Ref.  21).  In  particular,  as 
t  ->  0  (e.g.,  by  going  to  the  hard-sphere  limit  a  -*•  «,  F(2-rr (oj-cj )/av, )  -*•  1, 

w  U  •  I 

and  (2.35)  reduces  to  (2.15)1.  On  the  other  hand,  in  the  case  that  |u-w  |t 

v  W 

is  large,  specifically,  when 


27tJo)-{J 


av. 


»  1 


(2.37) 


(the  case  which,  as  will  be  seen  below,  obtains  in  the  far-wing  experiments 
in  C02),  the  "form-factor"  F(2tt(caj-oj0 )/av^ )  causes  a  sharp  reduction  in  the 
spectral  intensity,  I (to) ,  [relative  to  the  impact-theory  expression  (2.15)1. 
In  particular,  one  has,  to  a  sufficiently  good  approximation, 


/ 1  |  24tt2  -2tt  |  u)-o>0 1  /olv  ^ 


I(«)  * 


2ttt  (oj-u>0)‘ 


2  2 
a  v  ■ 


(2.38 


which  begins  to  look  like  the  empirical  formula  suggested  by  Benedict,  et  al.1 
However,  in  order  to  obtain  an  expression  which  can  be  compared  with  ex¬ 
periment,  the  rather  complicated  averaging  process  involving  the  consideration 

* 

of  all  types  of  impacts  (according  to  their  initial  trajectories)  as  well  as 
a  thermal  average  over  all  the  velocities  (those  of  translation  of  the  col¬ 
liding  partners  as  well  as  that  associated  with  the  rotatory  motion  of  the 
radiating  molecule).  The  systematic  treatment  of  the  total  averaging  process 
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will  be  deferred  for  the  future.  However,  certain  preliminary  studies,  in 
particular  those  concerned  with  the  thermal  averaging  of  the  relative  velocity, 
v^,  [appearing  In  the  exponential  factor  of  (2.38)1  will  now  be  presented.23 

The  basis  of  the  discussion  Is  the  fact  that  distribution  of  relative 
velocities,  v^,  is  given  by  the  standard  Loltzmann  relation 

l  -v^^Vy2 

P(vJ  =  7=e  1  1  ,  (2.39) 

1  VtirVj 

o  24 

where  v-j/  Is  the  mean  square  relative  velocity. 

The  normalization  coefficient  In  (2.39)  is  chosen  so  as  to  include 
both  signs  of  v^.  Since  we  are  momentarily  concerned  with  the  collision  of 
the  perturber  with  a  particular  "surface-element"  of  the  molecule  (eventually 
to  be  averaged)  we  consider  only  positive  v^  and  therefore  have 

2  -V^/SVy2 

p(v  )  =  -2  e  1  T  (2.39') 

1  ^7 

Inserting  (2.39')  Into  (2.38)  one  has 


I  (») 


1 

2ttt(w-u)q)2 


-2it|w-u0|/av 


(2.40) 
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where  <••••>  co] -|  Indicates  that  the  average  over  the  geometric  variables 
defining  collisions  has  yet  to  be  taken.  Looking  now  at  the  Integral  within 
the  braces  of  (2.40),  one  notes  that,  under  the  condition  specified  by  (2.37), 
the  exponential  factors  in  the  Integrand,  which  may  collectively  be  written 
as  e"F^  where  F(v)  =  2tt j o»-o)0 1 /av  +  v2/2vT2,  are  sharply  peaked  about  the 
value  v$  at  which  F(v)  becomes  a  minimum.  This  circumstance  suggests  the 
use  of  the  method  of  steepest  descents,  according  to  which  F(v)  is  approximated 
by  the  first  two  nonvanishing  terms  of  Its  Taylor  expansion  about  v$,  i.e., 


F(v)  -  F(vs)  +  \  F"(v$)  (v-v$)2 


(the  term  linear  in  (v-v$)  vanishing  since  v$  is  determined  by  the  condition 

that  F'(v  )  =  0);  moreover,  the  algebraic  factors  in  the  integrand,  such  as 

25 

Au>.j  and  v^,  are  approximated  by  their  values  at  the  saddle  point,  Au>$  and 
v$).  Omitting  the  details  of  the  calculation,  one  finds  that 


(a)  vs 

(b)  F(v$) 


3 

1 


-2/3 


(c)  F"(vs)  =  3/vT2 


so  that 


K«) 


(2.41) 


♦ 
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It  Is  now  of  Interest  to  note  that  the  ratio  l^l/Vg,  which  appears 
In  (2.41),  Is  In  fact  uniquely  determined  by  the  atomic  masses  of  the  col¬ 
lision  partners,  the  distance  R^q,  and  the  geometries  of  the  collision. 

Namely,  from  footnote  (20),  one  has,  for  any  collision, 

Au)  =  — y—  n  x  |  Fdt  a  ~y~  (n,  x  n^J  j | F | dt  ,  (2.42) 

where  ii  and  are  unit  vectors  along  the  momentary  molecular  axis  and  the 
llne-between-centers  of  the  colliding  atoms.  On  the  other  hand,  j[F|dt  Is 
directly  related  to  the  relative  velocity,  v;  In  fact 


(2.43) 
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where  M  Is  the  appropriate  reduced  mass  for  the  relative  motion  In  the 

2 

collision.  Its  calculation  must  be  equivalent  to  that  of  vT  ,  by  virtue  of 
the  fact  that  (2.39)  must  also  be  expressible  as 


1  _-Mvz/2kT  .  ,  _  ..  2  .  LT/U 

- — »  e  ;  i.e.,  vT  *  kt/m  . 

(2-vt2) 


An  explicit  expression  for  M  Is  given  In  Appendix  II  (cf.  text  below). 


Introducing,  then,  (2.42)  and  (2.43)  Into  (2.41)  one  has 
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where  use  has  been  made  of  the  fact  that  'co,  ■  2l&"o  <"o  *  mass  of  oxygen 
atom);  moreover,  from  Appendix  II  (Eqs.  2.24)  and  (2.25)],  one  has 


(kT/M) 


(2.45) 


M  “ 


1.1.. 

5  ^  ^ 


1  sin20 


(2.46) 


where  M^,  MCq  ,  and  Mg  are  respectively  the  masses  of  the  noble-gas  perturber- 
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atom,  the  C02  molecule,  and  an  oxygen  atom  and  6  the  angle  between  the 
molecular  axis  and  the  llne-of-centers  of  colliding  atoms  (sine  *  |n  x  nj). 
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3.  COMPARISON  WITH  EXPERIMENT  AND  DISCUSSION 

Let  us  begin  by  computing  an  upper  bound  for  the  magnitude  of  (2.44) 

at  a  frequency  u  *  2600  cr'*,  such  that  |w-w0l  *  250  cm'*.  This  region 

corresponds  to  the  maximum  frequency  attained  In  the  cited  experiments 

(cf.  Ref.  1,  Fig.  3,  and  Ref.  2,  Fig.  14)  on  the  fundamental  v3  band  of  C02- 

The  motivation  for  such  a  computation  Is  our  anticipation  that  the  result 

will  fall  substantially  below  the  experimental  values.  If  this  be  the  case, 

we  shall  regard  the  situation  as  promising,  since  we  have  not  yet  Invoked 

attractive  Interactions;  according  to  preliminary  calculations  these  should 

28 

give  rise  to  an  Increase  In  absorption. 

It  will  be  noted  that,  apart  from  the  factor  2  occurring  In  its  first 
denominator,  (2.44)  Is  the  product  of  a  Lorentzlan  line  shape,  multiplied  by 

the  quantity  ^ . }coir  Although,  as  has  been  remarked  In  Section  2 

Icf.  text  subsequent  to  Eq.  (2.16)1.  the  factor  1/t,  occurring  In  (2.44), 

Is  expected  to  be  substantially  smaller  than  the  corresponding  quantities 
occurring  In  the  standard  Lorentz-type  formulae  (which  apply  to  the  near 
vicinities  of  Individual  rotational  components),  we  shall  ignore  this 
difference  for  the  time  being.  The  main  discrepancy  between  the  Lorentz- 
type  extrapolations  used  In  the  cited  experiments  (cf.  Ref.  1,  Fig.  2)  and 
the  experimental  curves  will  therefore  be  sought  In  the  factor, 
of  (2.44),  which  we  now  proceed  to  estimate. 

For  numerical  values  of  the  various  quantities  occurring  In  (2.44) 
(and  In  the  subsidiary  expressions  (2.45)  and  (2.46)  for  Vj  and  M),  we  take 
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a  * 

4  x  108  cm*1 

(3.1) 

(ii-ii)  * 

0 

2ttc (250)  -  (1.9  x 

1011)  (250)  sec*1 

(3.2) 

kT  « 

300  (1.38)  x  10‘16 

-  4.15  x  10’14ergs 

(3.3) 

ma  * 

40  Mh 

(3.4) 

"co2  * 

44  M« 

(3.5) 

"0  8 

16  "h 

(3.6) 

Rco  8 

10*8  cm*1 

(3.7) 

where  ^  *  1.67  x 

-24 

10  gm  is  the  mass 

of  the  hydrogen  atom.  Utilizing 

(3.3,  4,  5,  6),  one  finds  (upon  setting  sine  «  1;  It  is  shown  below  that 

this  choice  maximizes  <(*•••  -)coi -j ) 

M  - 

12.7  H„ 

(3.8) 

vT  * 

1.59  x  107(12 .7  )** 

a 

«  4.5  x  10  cm/ sec  . 

(3.9) 

From  these  values,  one  has 


so  that 


9.75 


(  •  *>con 


(2.45)  (.63)  (4/3)*1  e*9*75 
1.78  e*9* 75  •  (1.78)  (.58  x  10*4) 


(3.10) 

(3.11) 


1.03  x  10*4  .  (3.12) 
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Let  us  now  note  that  setting  sine  s  1  minimizes  M,  thereby  maximizing 
vT;  the  result  Is  the  maximization  of  the  exponential  factor  In  (2.44).  The 
remaining  sine-dependence,  contained  In  the  pre-exponential  factor,  may  be 
written  as  the  square  of 


_ sine _ 

W  V"Co2*  0/2)s1nze 


sine 

.764  ♦  .5  s1nZe 


(3.13) 


which  Is  easily  shown  to  be  a  monotonlcally-increasing  function  of  sine  in 

the  range  0  i  sine  £  1.  Thus  ^ . )coll  maximized  by  taking  sine  equal 

to  unity,  q.e.d. 

To  proceed  further,  we  may  take  the  liberty  of  using  for  1/t  a  value 
equal  to  that  quoted  by  Benedict  et  al  (cf.  reference  1,  page  536)  for 
CO2-  C02  collisions,  namely  one  ccrrespondl ng  to  a  width  parameter  y  *  1/2ttct, 
of  .094  cm” *  (at  one  atmosphere,  STP;  this  condition  Is  used  herein  and  in 
what  follows  for  all  density-dependent  quantities.  Such  a  value  for  1/t  is 
quite  maximal,  first  of  all  because  of  the  remarks  in  the  second  paragraph 
of  this  section,  and  also  because  the  y  for  self-broadening  Is  known  to  exceed 
that  for  foreign-gas  broadening  by  a  factor  -v  1.5  (cf.  reference  1,  page  536; 
also  reference  2,  Table  III). 

The  absorption  coefficient,  k(u),  is  given  by  the  expression 

k(u)  •  Sl(u>)  ,  (3.14) 
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where  S  Is  the  total  line  strength;  Its  value  for  the  fundamental  band 
of  CO2  (as  given  by  reference  1,  bottom  of  page  535)  Is 

S  *  2700  cm”*  atm”*. 

Substituting  (2.44)  Into  (3.14),  one  has2' 

k(u>)  *  - - - 7  ('*•••) coll 

2ttt(u>-u0) 

1>03x10-4 

2*(250r 

=  .65  x  10“2  (1.02  x  10”4)  =  .67  x  10”7cm”*  atm”2  , 

which  Is  lower  than  the  value  given  by  Burch's  measurement  for  argon  (cf. 

reference  2,  Fig.  10)  by  a  factor20  %  15.  In  view  of  all  the  maximal  estimates 

invv/.ved  in  the  above  calculation,  it  seems  apparent  that  a  new  physical 
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ingredient  has  to  be  introduced  into  the  theory. 

However,  as  remarked  in  the  first  paragraph,  such  an  Ingredient  already 
Impends:  namely,  the  existence  of  attractive  Interactions,  especially  the 
rapidly  varying  ones,  such  as  the  dispersion  (van-der-Waals)  interactions 
*  r”6,  which  attain  considerable  magnitudes  at  the  distance  of  closest  ap¬ 
proach.  Preliminary  exploration  Indicates  that  these  interactions  can  be 
incorporated  into  the  mathematical  apparatus  of  the  theory.  It  would  there¬ 
fore  be  highly  desirable  (and  is  in  fact  intended)  to  pursue  this  direction 
in  the  near  future. 


4. 


FUTURE  DEVELOPMENTS  AND  PROSPECTS 


The  concluding  remarks  of  the  preceding  section  have  brought  us  to  the 
stage  where  a  discussion  of  further  developments  and  overall  prospects  is  now 
in  order. 

A.  Further  Developments  of  the  Present  Theory 

1.  As  pointed  out  in  the  last  section,  the  most  urgent  requirement 
is  the  introduction  c(  attractive  interactions;  it  is  expected  (cf.  footnote  31) 
the  effect  of  these  would  be  to  shorten  xc  and,  hence,  to  increase  the  magni¬ 
tude,  while  decreasing  the  slope,  of  I (cu)  in  the  far  wing. 

2.  A  point  of  more  theoretical  (interpretative) interest  is  the 

feature,  noted  in  the  excursus  following  Eq.  (2.20)  of  Section  2,  that  the 

frequency-dependence  of  our  impact  result  for  I((d)(<*  1/(u)-<jj0)4) differs  from 

the  Lorentz-type  vari&cion  (<*  l/(w-w0)2},  usually  exhibited  by  impact-type 
13 

theor’es.  A  tentative  explanation  was  suggested  in  the  text  subsequent  to 

Eq.  (2.20).  It  remains  to  verify  this  suggestion  by  a  detailed  analysis; 

from  it  we  may  hope  to  achieve  a  deeper  physical  understanding  of  the  standard 

13 

Lorentz-type  theories,  based  on  the  cited  references, 

3.  The  preliminary  estimates  of  the  smallness  of  Ivl-b(w)  [which 
denotes  the  contribution  of  vibrational  perturbations  to  I(w)],  mentioned  at 
the  end  of  the  Introduction,  should  be  reinforced  by  a  detailed,  systematic 
investigation.  In  particular,  it  will  be  desirable  to  alter  the  treatment  of 
reference  4  (which  is  for  the  case  of  a  diatomic  molecule)  to  apply  to  the 
somewhat  more  complicated  situation  of  a  linear  tri atomic  molecule,  such  as  CO 

4.  Since  the  maximum  values  of ti (w-w  ) ,  encountered  experimentally, 
are  not  really  small  compared  to  kT,  quantum  modifications  of  our  treatment 
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ro 


should  be  considered  In  this  connection,  it  may  be  remarked  that  the  treat¬ 
ment  of  given  in  reference  4  is  quantum-mechanical,  permitting  a 

comparative  study  ffor  arbitrary  fi(u>-u)0)]  with  the  classical  limit  [attained 
when  1Ku)-cu0)  <<  kT]  It  is  proposed  to  carry  out  such  a  study;  if  it  shows 
that  quantum  modifications  are  important  under  experimental  conditions,  an 
attempt  to  introduce  them  into  the  present  theory  should  he  made. 

It  must  however  be  stated  that,  even  with  all  the  above-proposed 
improvements  -  including,  of  course,  the  proper  evaluation  of  l/i  and 
‘ ‘Xoll  -  1,e  *  a  Pr°Per  averaging  over  the  geometries  of  collisions 
(in  place  of  the  upper-bound  estimates  given  in  Section  3),  the  theory  will 
still  contain  residual  uncertainties,  associated  with  those  of  intermolecular 
force  laws;  such  uncertainties  are  augmented  in  I(u>),  due  to  the  exponential 
dependence  of  this  quantity  on  interaction  parameters  (as  shown,  e.g.,  by 
(2.41);  preliminary  indications  are  that  the  introduction  of  attractive  inter¬ 
actions  will  not  affect  this  feature).  Looking  at  the  situation  positively, 
however,  one  may  utilize  the  present  theory  [augmented,  of  course,  by  the 
proposed  extensions  and  improvements  (especially  the  introduction  of  attractive 

forces))  as  another  method  for  the  determination  of  intermolecular  force  laws  - 
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to  be  used  in  conjunction  with  various  other  methods. 

From  the  practical,  or  engineering  point  of  view,  a  proper  deter¬ 
mination  of  the  interaction  parameters  contained  in  I (to)  will  permit  extra¬ 
polation  to  frequency  regions  which  have  not  as  yet  been  observed  experimentally. 
One  would  thus  have  available  a  useful  tool  for  the  study  of  absorption  in  the 
so-called  atmospheric  "windows,"  in  which  absorption  is  usually  dominated  by 
the  far-wings  of  strong  bands 
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B.  There  remains  one  further  subject  for  future  research,  which  has 

not  been  discussed  as  yet.  Namely,  as  was  noted  in  the  introduction  (Section  1), 
2 

Burch  et  al  found  that,  in  the  high-frequency  far  wing  of  the  self-broadened 
3v3  band  of  CC^,  at  7000  cm'1,  the  absorption  coefficient  k(w),  although  sub- 
Lorentzian,  is  much  higher  than  would  be  inferred  from  the  comparable  situation 
in  the  case  of  the  fundamental  band  in  the  vicinity  of  2400  cm"1.  In  what 
follows,  certain  speculations,  which  in  our  opinion  are  relevant  for  the  ex¬ 
planation  of  this  rather  unique  phenomenon,  will  be  outlined.  Briefly,  we 
consider  the  two-stage  process  described  by  the  questions, 


CO 

CO 


+  tia)  -*■  CO 


(3) 


+  CO 


<°>  *  co<2> 


In  these  equations  the  superscript  attached  to  CO,,  indicates  the  vibrational 
quantum  number  of  the  v3  branch  (those  of  the  other  branches  being  zero 
throughout).  The  first  stage  occurs  via  electromagnetic  interaction;  the 
second  may  occur  via  dipole-dipole  coupling  (which  is  strongly  allowed  by 
harmonic-oscillator  selection  rules);  also,  the  adjective  "two-stage"  (used 
above)  implies  that  intermediate  energy  need  not  be  conserved. 

Now,  in  considering  the  physical  significance  of  the  above  mechanism, 
we  note  the  following  important  fact.  Namely,  the  exponential -type  fall-off 
characteristic  of  the  present  theory  depends  crucially  upon  the  discrepancy 
between  the  energy  of  the  absorbed  light  quantum  and  the  vibrational  energy 
of  the  absorbing  system  (as  is  clearly  put  in  evidence  by  the  occurrence  of 
the  quantity  (u)-oj0 )  in  the  argument  of  the  exponential  factor  in  the  formulae 
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of  Section  2)  With  this  fundamental  point  in  mind,  let  is  now  note  that, 
due  to  small  anharmonicity,  the  vibrational  energy  of  the  final  state 
CO ^  *  CO^  turns  out  to  be  larger  than  that  of  intermediate  state  CO^; 
in  fact,  from  data  in  Herzberg's  book,  one  finds  that 

1^20  q  ”  ^30  ^  50  cm 
so  that  in  the  high  frequency  far  wing 

1i[w-  (a>20  +  u10)]  <  ^(«“«30)  » 

i,e,,  the  two-stage  reaction  permits  the  realization  of  a  smaller  discrepancy 
between  the  energy  of  the  absorbed  quantum  and  that  of  the  final  vibrational 
state.  From  this  result  we  are  inclined  to  infer  that  the  onset  of  the  ex¬ 
ponential  fall -off  would  occur  at  higher  values  of  u>-u>0  than  would  be  ex¬ 
pected  in  the  absence  of  this  mechanism.  As  may  be  seen  from  Fig,  6, 
reference  2,  this  spectral  behavior  corresponds  to  the  experimental  observations. 

The  above  remarks  constitute  the  bare  framework  of  an  idea  which 
we  propose  to  use  as  a  basis  for  the  development  of  a  theory  which  we  hope 
will  provide  an  explanation  of  the  anomalous  absorption  observed  by  Burch 
for  th .  rase.  In  this  connection,  we  have  been  encouraged  by  recent  infor¬ 
mation  from  Burch  (private  conversation)  that  similar  behavior  has  been 
observed  for  the  high  frequency,  far  wing  of  the  self-broadened  3v^  line  of 
N20,  Examination  of  the  relevant  anharmonicity  parameters  in  Herzberg's 
book  indicates  that  here  also  the  above  described,  two-stage  process  is 
energetically  favored 
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APPENDIX  I 


4 

4  » 


In  this  appendix  the  integral 

.+00 


J.  - 


o  aVjt  -i  (w-ai  )t 
sech2  -i-a  0  dt 


(1.1) 


will  be  evaluated.  To  this  end  it  is  first  convenient  to  introduce  a  new 
variable  of  integration 


av.jt/2 


and  rewrite  (1.1)  as 


J.  =  1 

i  2 


+00 


o  -2i(a)-u)  )u/av. 
sech^C  e  0  1 


(1.2) 


One  then  considers  a  clockwise  integration  over  the  rectangular  contour  whose 

corners  are  the  points  -<*>,  +«,  +®-iir,  -«-iir.  The  integration  over  the  top 

side  (-»  +  +«)  gives  (1.2).  For  the  integration  over  the  bottom  side 

2 

(-Itt  -*■  -oo+itt),  one  notes  that  (a)  the  sech  ^  factor  remains  unchanged,  (b) 

-2tt(u)-u)  )/av . 

the  exponential  factor  gets  multiplied  by  a  constant  e  1  and  the 

-  2  tt(  u)-u) l  )/av. 

sign  is  negative;  the  result  is  then  -e  .  Since  the  contributions 

2 

of  the  vertical  sides  are  zero  (by  virtue  of  the  sech  £  factor), one  has 


j 


h 


1-  e 


1 

7 


-2ir(w-u)0)/av^  *1  ’ 

o  -2tt(u-u  )£/av, 
secfre  e  0 


(I  3) 


d£ 
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where  Cp  denotes  a  counter-clockwise  Integration  around  a  small  circle  the 
center  of  which  Is  the  single  pole  of  the  Integrand  (within  the  original 
rectangle)  at  the  point  £  =  -tt1/2.  Expanding  sech  £  In  the  neighborhood  of 
£  =  --n-i/2,  one  has  (with  £  =  -  +  z) 


,  -tt(o)-u)  )/av,  t  -2i(w--J0)z/avi 

ie  0  J- - 7 - 


2ir(a)-a)o)  -7r(u)-u0)/a ) 


a  V- 


(1.4) 


which,  together  with  (1.3),  yields 


J 


av, 


1-e 


■2ff(w-uft)/av 


1 


(1.5) 


APPENDIX  II 


The  subject  matter  of  this  appendix  is  the  derivation  of  an  expression 
for  the  probability  distribution  of  the  relative  velocity  of  two  colliding 
atoms,  one  of  which  is  a  constituent  of  a  triatomic  linear  molecule.  For  the 
sake  of  definiteness  we  consider  the  text  case  of  a  noble-gas  atom,  e.g., 
argon,  colliding  with  one  of  the  oxygen  atoms  of  C02.  One  task  is  then  to 
derive  (2.39),  i.e.,33 


P(v) 


(II. 1) 


where  v  is  the  component  of  relative  velocity  of  colliding  noble  gas  and 

2 

xoygen  atom,  together  with  an  explicit  expression  for  vT;  the  calculation 
will  also  yield  [cf. text  subsequent  to  Eq.  (2.43)]  an  expression  for  the 
reduced  mass,  M. 

The  starting  point  of  the  calculation  is  an  expression  for  the 
probability  distribution  of  the  various  elementary  velocities  involved  in 
the  problem,  to  wit 


Pl»A.  V  *r>  d W2ir 


-*Xa 

- 37?  ' 


l2”vTA 


d\  -^Tg 

tV  ' 


lK,» 


d^v  -v^/>v t 

— r  r/w  Tr 

- 2 —  e 


2ttv 


Tr 


(II- 2) 
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In  this  expression,  v^,  v^  and  v^.  respectively  denote  the  vector  velocities 

of  the  noble-gas  atom,  the  center  of  gravity  of  the  COg  molecule,  and  the 

3  3 

rotational  velocity  of  said  molecule.  The  notations  d  v^  and  d  v^  refer 

to  the  three-dimensional  (cartesian)  volume-elements  of  v^  and  v^;  the 
2 

notation  d  v^.  takes  account  of  the  fact  that,  by  virtue  of  the  condition 

v/n  =  0  ,  (II. 3) 

~T  — 

there  are  only  two  independent  components  of  rotational  velocity,  both 
perpendicular  to  n 

On  the  r.h.s.  of  (I I. 2),  the  constants  VjA,  vTg,  and  vTr  are  given  by 
the  formulae 


VTA  * 

kT'HA 

(II. 4a) 

vTg  s 

kT/Mg  = 

kT'"co2 

(II. 4b) 

vTr  = 

-\i~ 

II 

kT 

(II.4c) 

In  the  case  of  the  first  two  square  brackets,  the  insertion  of  (II  4a) 
and  (II  4b)  yield  the  manifestly  correct  Boltzmann  distributions  (the  factors 
|£tvjJ  and  {2^]  taking  care  of  the  normalizations). 

The  third  square  bracket,  describing  the  distribution  of  the  rotational 
velocity,  requires  some  comment.  First  of  all,  we  note  the  obvious  feature 
of  its  difference  in  form  from  the  other  two  brackets;  this  difference  arises 
solely  from  the  above-noted  two-dimensional  character  of  the  motion.  In  order 
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to  clarify  the  bracket  further,  we  note  that  It  may  be  written  in  the  form 


P(u))d  CO 


2  IT. 


*RC0‘‘)2/2v 


2 

Tr 


(w  =  angular  velocity) 

which,  with  the  use  (II  4c),  becomes 


d2a  e-Iu2,'2kT 


with 


2 

Tr 


=  vTr^RC0 


(II- 5) 


«  kT/I  ,  (II. 6) 

2 

and  where  the  notation  d  u  takes  note  of  the  fact  that  w  is  also  a  two- 
dimensional  vector,  perpendicular  to  the  molecular  axis,  n. 

Equation  (II. 5)  already  begins  to  resemble  the  standard  expression; 

-IwJ/2k_ 

in  particular  with  respect  to  the  Boltzmann  factor,  e  In  order  to 

achieve  an  even  more  familiar  form,  let  us  go  over  to  a  two-dimensional 

polar  coordinate  system,  in  which  u  is  expressed  in  terms  of  its  magnitude, 

and  angle,  $>,  with  respect  to  an  axis  perpendicular  to  in,  but  otherwise 

2 

arbitrary;  then,  upon  using  the  formula  d  w  *  wdua$  and  integrating  over 
one  obtains  the  standard  angular-velocity  distribution  function 

P(*)d*  *  e-lu2/2kT  t  (11.7) 
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or.  In  terms  of  angular  momentum,  J  =  Iw 


P(J)dJ 


JdJ 

TFT 


e-J2/2IkT 


(II. 8) 


It  should  finally  be  remarked  that  the  third  bracket  of  the  r.h.s.  of 
(II. 2)  may  be  obtained  from  tt.c  fundamental  prescription  of  statistical 
mechanics,  which,  for  the  case  of  rotational  motion,  reads 

1  -vkT 

P(e,*;pe,p^)d*dedp^dpe  *  -  e  d^dedp^  .  (ILj 


where 

Hr  *  2t[p8  *  P»/S,n2e 

Is  the  rotational  Hamiltonian,  d$d6dpjto„  is  a  differential  volume  element 
in  the  four-dimensional  phase  defined  by  the  rotational  coordinates  6,  <J>,  and 
associated  canonical  momenta,  pQ  and  p^,  and  z  is  the  partition  function 
(which  will  not  be  written  down  explicitly  here,  since  from  our  point  of  view 
it  is  merely  a  normalizing  coefficient).  Briefly,  the  procedure  is  transform 
from  the  canonical  momenta  to  the  associated  velocities,  via  the  formulae 


Vg  =  Rcq0  =  pe 


(II  11a) 


\  5  RC0{sin6)*  -  TsT^P*  • 


(II. lib) 


thereby  obtaining 
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P(e.$;pft.pjded*dp  dp. 


,2  -U/R*,)[.^]/2kT 


(II  12) 


zr: 


CO 
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Note,  now  that  (11.12)  is  equal  to  the  typical  spatial  volume  element 
dfln  =  sm0d0d$  (defining  the  direction  of  n)  multiplied  by  a  factor  which, 
in  view  of  (II  4c),  may  be  written  as 


P(vQ,^)dv0dv^  =  (const)  e 


21  2 
Tr 


dvedv* 


(II  13) 


Since  the  direction  of  v0  and  v^  are  both  perpendicular  to  n  (as  well  as  to 
each  other),  the  r.h.s.  of  (11.13)  is  in  fact  equal  (apart  from  a  normalization 
constant)  to  the  third  bracket  of  (II. 2),  q.e.d. 

With  these  preliminary  remarks  on  the  validity  of  (II. 2)  out  of  the 
way,  we  now  proceed  to  a  second  basic  formula,  namely,  the  relationship 
between  v  and  the  vector  velocities  occurring  in  (I I. 2).  Remembering  that  v 
denotes  the  component  of  relative  velocity  of  colliding  atoms  (in  our  example 
the  noble-gas  perturber  h  one  of  the  oxygen  atoms  of  C02)  parallel  tc  the 
line-of-centers  of  said  c  jms  (i.e.,  to  r^),  we  have 

v  *  * 

which,  upon  going  ever  to  a  cartesian  system  whose  z-axis  is  parallel  to  the 
molecular  axis,  n,  may  be  written  as 

'  •  *  "cy'WV 

*"«(,ArV’  (U  111 

one  notes  the  absence  of  the  term  linear  in  vri  since,  with  n  along  the  z-j*:s 
(II. 3)  becomes  vri  *  0. 
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In  proceeding  further,  it  is  expedient  to  go  over  to  a  more  general 


notation,  In  which  (II. 2)  Is  written  <n  the  form 


P(  —  Vj  —  )dvj. . .  dvj . . . 


8  f  ^ 

=  n  - - 

1  =  1  (2rrv|1 


2  2  • 

■Wn 


(11.15) 


where  1  refers  to  the  cartesian  components  of  the  various  velocities.  For 
the  sake  of  definiteness,  1  -  1,2,3  denote  the  x,y,z  components  of  v^, 
i  =  3,4,5  the  corresponding  components  of_v^,  and  i  =  7,8  tne  x,y  components 
of  jvr  Inspection  of  relations  (II.4a,b,c)  then  show  that 


VT1  =  VT2  '  VT3  s  kT/KA 


VT6  *  kY/MC0, 


kT/2M,. 


In  this  new  notation  (11.14)  may  be  written  as 


(II. 16a) 


(II.16b) 


(II  16c) 


v  =  1 


r<  . 


where 


al  “  *4  *  a7  *  "c* 


a2  *  a5  “  *8  “  nc 


a3  *  a6  *  ncz  * 


(II  17) 


(II.18a) 


(II  18b) 


(11.18c) 


r 


P. 


k 


I 


The  normalized  expression  for  P(v)  is  now  written  as 


P(v)  = 


dv, 


dv8  P( - vi - )6(v-X  a.v.) 

i 


(11,19) 


Equation  (11,19)  expresses  the  fact  that  P(v)  is  qotten  by  integration  over 
all  values  of  v.  which  are  subject  to  the  constraint  imposed  by  (11.17) 

(and  incorporated  in  (11.19)  by  the  delta-function).  Expressing  the  delta- 
function  in  its  Dirichlet  representation 


M*>  ■  £ 


+OC 


Jxt  ,4  + 
e  dt 


and  utilizing  (11.15),  one  then  has 


P(v)  = 


1 


+oo 

JVt 

dt  e  it. 


+00 


dv. 


-V2v?i 


’taivi 


(2iTVji 


(II  20) 


The  evaluation  of  the  veloc-tv  integrals  is  standard;  one  obtains  for  each 

‘t’vTiai 

square  bracket  in  (11.20)  a  factor  e  so  that  (11.20)  becomes 


P(v)  = 


_1_ 

27! 


+”  -tV/2 

ivt  x  VT^ 
e  e  dt 


(11-21) 


where 


v 


2 

T 


l 


(II  22) 


The  t-integration  is  now  done,  yielding. 


i  -v  /2vt 

p(v)  =  — irr  e  (11.23) 

which  is  identical  to  (II. 1),  q.e.d.  Moreover,  the  required  explicit  expression 
for  vT  is  given  by  (11.22),  in  conjunction  with  formulae  (II.18a,b,c)  and 
(II,16a,b,c).  Combining  these  relationships  one  has,  without  further  ado 


v 


2 

T 


1  .  2 

^Sln  9 


(11.24) 


where,  with  e  the  angle  between  nc  and  the  z  axis  -  that  is,  between  _n^  and 
the  molecular  axis  n,  the  obvious  relation 


.  ? 

Sin'G  = 


"cx  + 


cy 


=  (n  x  n CY 


has  been  used. 

We  may  now  note  briefly  that  (11.24)  used  in  conuunction  with  the  text 
relationship 

v^  =  kT/M 

yields  for  the  reduced  mass,  M,  the  expression 
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16.  The  time-dependence  of  this  relative  velocity  will  shortly  be  shown  to 

be  intimately  linked  up  with  (in  fact  essentially  identical  to)  that  of 

ul  in  a  collision  interval. 

-P 

17.  The  r.h.s.  of  (2.25)  is  obtained  from  the  more  basic  relation 

v  -  (2/M)^  (E  -  V  e'ax j*5 

(where  M  is  the  appropriate  reduced  mass)  by  factoring  from  the  square 
root  the  initial  "normal"  energy,  E  =  |  Mv^,  and  introducing  xc  =  ^log  ^  . 

18.  In  the  case  of  a  noble  gas  of  intermediate  mass,  such  as  argon,  one  may 
take  vrot/v  v  1/vf;  the  smallness  parameter  is  then  ^  1/5-6. 


19. 


Incidentally,  pre 
are  quadratic  in 


iminary  considerations  indicate  that  said  corrections 

vrot  1  \ 
v  2R'C0  J  * 


20,  Within  the  framework  of  (2,24)  and  (2.27),  the  overall  frequency  change, 
Auj^  ,  is  simply  equal  to  (Rqq/I)^  x  j  Fdt;  explicit  values  may  always 
be  obtained  by  going  to  the  impact  limit  (a  -*  «>)  and  using  standard 
methods  of  classical  mechanics  (as  is,  in  fact,  done  in  the  latter  part 
of  this  section) 


21.  It  will  be  noted  that  the  factor  (t1)  sech2[  av. (t-t^ )/2]  reduces  to 
6 (t~t  1 )  in  the  limn  a  -*■  *  ,  i  .e. ,  in  the  limit  in  which  the  repulsive 
potential  becomes  that  between  two  hard  spheres  For  such  a  potential 
collisions  are  clearly  pure  impacts.  In  the  case  of  finite  a,  it  is 
also  clear  that  7C  £  2/av. 


* 
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22.  This  statement  is  based  on  preliminary  studies  which  also  indicate 

that  the  contribution  of  ( . )..  is  quadratic  in  the  new  parameter 

A  systematic  study  is  reserved  for  the  future, 

23.  Since  vi  appears  in  the  rapidly  varying  exponential  factor  of  (2.38), 

we  believe  that  the  thermal  average  is  the  most  important.  In  particular, 
as  will  be  seen  below,  it  leads  to  an  exponential  factor  whose  argument 
is  fractional  in  (o)-cj0),  rather  than  linear  [as  is  superficially  ex¬ 
hibited  by  (2,38)1. 

24.  More  generally,  it  is  being  asserted  that  the  distribution  of  the  pro¬ 
jection  of  an  arbitrary  sum  of  velocities  along  an  arbitrary  direction 
obeys  (2. 39),  with  a  suitably  chosen  v-j-.  This  assertion  is  easily  proved 
for  the  case  when  all  the  component  velocities  refer  to  purely  linear 
motion.  Its  proof  for  the  case  when  one  of  the  components  is  rotational 
is  somewhat  more  subtle,  and  will  be  given  in  Appendix  II.  It  is  also 
subject  to  the  limitation  of  (2  27);  only  within  that  limitation  is  it 
possible  to  regard  the  rotational  motion  as  linear, 

25  It  is  here  presumed  that  the  rotational  velocity  is  also  involved  m 
the  thermal  average 

26-  All  other  masses  have  appropriate  subscripts  attached  to  them,  so  no 
notational  problems  arise. 
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27, 


Under  the  assumption  [incorporated  by  Eq.  (2,23)]  that  the  Interaction 
between  the  perturber  and  the  molecule  is  a  sum  of  spherical  interaction 
between  the  perturber  and  the  individual  atomic  constituents,  and  with 
the  condition  expressed  by  (2.24),  it  is  clear  that  only  perturber- 
oxygen  collisions  contribute  to  nonvanishing  Au^'s. 

28.  A  similar  effect  is  known  to  exist  for  the  very  analogous  case  of  vibra¬ 
tional  deactivation  (cf.  K.f,  Herzfeld  and  T.A.  Litovitz,  Absorption  and 
Dispersion  of  Ultrasonic  Waves,  Acad.  Press  Inc,,  New  York,  1959, 

pp  278-285 

29.  If  one  expresses  l/i  ir  wave  numbers  (which  is  equivalent  to  replacing 

it  by  y)  and  expresses  frequencies  also  in  wave  numbers,  k(oo)  comes  out 
-1  -2 

in  units  of  cm  atm  (STP),  as  is  to  be  expected  since  it  is  proportional 
to  the  product  of  densities  of  absorbing  and  perturbing  gases, 

30.  Alternatively,  we  can  recalculate  the  theoretical  k(co)  by  taking  w-u)  = 

210  cm’*  (corresponding  to  ui  •-  2560  cm"*,  at  which  an  actual  measurement 
was  made);  it  turns  out  that  gets  increased  by  a  factor  of 

(2  86),  which  together  with  a  factor  (250/21 0)2  (which  accounts  for  the 
change  m  (w-*  )  ),  yields  a  theoretical  k(w)  of  2.7  x  10  cm  atm  , 

to  be  compared  with  Bunch's  value,  -  2  4  x  10  cm  atm’  Thus  the 

theoretical  kU)  is  still  lower  than  the  experimental  value  by  a  sub¬ 
stantial  factor  (^  10). 
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31.  It  may  also  be  remarked  that  the  average  theoretical  slope  of 

logg  ^ _ )coll  vs  ln  ran9e  '3etween  210  cm’1  and  250  cm’1 

i  ’ 

Is  k  .026  cm,  whereas  the  corresponding  experimental  slope  (as  closely 
as  we  are  able  to  estimate  it  from  Fig.  14  of  Ref.  2)  is  %  .014  cm. 
However,  in  this  region,  experimental  uncertainties  are  substantial; 
moreover,  (and  more  importantly)  the  theory  in  its  present  form  requires 
much  further  development  (and  refinement)  for  a  comparison  of  such 
details  as  slopes  to  be  of  any  real  significance. 

32.  For  a  review  of  these,  see  Advances  in  Chemical  Physics,  Vol ,  12: 
Intermolecular  Forces,  edited  by  J.O.  Hirschfelder  (Interscience 
Publishers,  1967) 

33.  The  subscript,  i,  in  (2.39),  is  dropped  for  convenience. 

34.  Note  that  in  relating  the  differentials  of  the  momenta  to  those  of 
tne  associated  velocities,  the  coordinates  <j>  and  6  are  held  constant 
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CHAPTER  II 
BEAM  OYNAMICS 
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I.  INTRODUCTION 

This  report  is  concerned  with  the  effect  of  target  motion  on  thermal 
self-defocusing.  In  a  previous  study^  we  have  produced  a  description  of  the 
self-defocusing  effect  applicable  only  to  the  case  of  a  stationary  target, 
that  is,  to  the  situation  in  which  the  laser  beam  heats  a  fixed  air  mass 
throughout  the  entire  pulse  period.  If  the  target  is  moving,  however,  the 
beam  must  be  rotated,  and  a  given  air  mass  is  heated  only  for  the  time 
required  for  the  beam  to  pass  through  it;  consequently,  there  must  be  less 
defocus ing. 

Jt  was  found  that  for  a  stationary  target  the  limitation  caused  by 

self-defocusing  is  on  the  total  amount  of  energy  that  can  be  delivered  to 

the  target.  However,  *or  a  moving  target  we  expect  the  limitation  to  be  on 

the  deliverable  power.  The  reason  there  is  no  limitation  on  the  total  enerqy 

is  simply  that  new  air  masses  are  constantly  being  used  for  the  beam  path, 

so  there  is  no  cumulative  heating.  By  analogy  to  the  stationary-target  case 

we  expect  a  limitation  on  the  energy  deliverable  through  any  given  air  mass, 

so  the  maximum  transmi ttable  power  depends  on  how  rapidly  the  beam  rotates. 

The  amount  of  time  a  particular  air  mass  is  heated  depends  on  its  distance 

from  the  laser,  and  is  smallest  at  the  target  distance  (its  value  there  is 

simply  the  beam  width  divided  by  the  target  velocity).  Therefore,  we  can 

obtain  a  simpij  overestimate  of  the  maximum  transmittable  power  from  the 

% 

energy  limit  calculated  in  the  stationary  target  case1,  namely,  the  latter 
energy  divided  by  the  heating  time  at  the  target. 

In  this  report  we  obtain  a  more  precise  estimate  of  the  maximum  trans- 
mittable  power  by  extending  the  stationary-target,  analysis.  The  extension  is 
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somewhat  crude,  but  it  at  least  accounts  for  the  dependence  of  the  heating 
time  of  any  air  mass  on  its  distance  from  the  laser.  The  crudity  lies  in  the 
way  that  we  treat  the  beam  rotation.  Instead  of  solving  equations  that 
properly  describe  a  beam  continuously  rotating  into  new  regions  of  air,  we 
use  the  equations  that  describe  a  stationary  beam;  we  simulate  the  rotation 
by  assigning  to  the  air  at  any  given  distance  from  the  laser  a  heatlnq  In¬ 
terval  estimated  from  the  time  it  would  take  a  rotating  beam  to  pass  through 
it 

The  stationary-beam  analysis  is  based  on  the  ray  optics  approximation, 
that  is,  diffraction  spreading  is  ignored.  In  addition,  it  is  assumed  that 
the  time  required  for  the  beam  to  pass  through  any  air  mass  is  longer  than 
the  time  required  for  sound  to  propagate  in  air  across  the  beam.  Further¬ 
more,  another  approximation  is  made,  which  significantly  simplifies  the 
mathematical  system.  An  integral  which  appears  in  the  equations,  namely  the 
time  integral  of  the  light  flux  at  any  point  in  space,  is  replaced  by  the 
simple  product  of  the  instantaneous  flux  and  the  time.  As  a  result  the  time 
becomes  a  trivial  parameter  instead  of  an  independent  variable,  and  the 
equations  become  soluble  for  a  certain  special  case  Now,  in  the  extension 
described  in  this  report,  what  we  have  done  to  simulate  the  beam  rotation  is 
to  replace  the  actual  time 'in  this  product  of  flux  and  time)  by  the  time 
interval  over  which  the  air  mass  at  the  point  in  question  has  been  illuminated. 
Although  this  time  interval  now  depends  on  one  of  the  independent  variables, 
namely ,  the  distance  ?*■©*  c;-e  laser,  we  find  that  tie  equations  are  still 
soluble  for  the  sane  special  case  as  before. 
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The  particular  manner  in  which  the  beam  spread*;  by  self-defocusing  has 
been  found  to  depend  on  the  specific  shape  of  its  initial  power  profile.  For 
example,  if  the  power  profile  is  uniform  in  the  center  and  falls  off  only  at 
the  edges,  then  initially  the  rays  in  the  uniform  interior  will  not  deflect 
at  all,  but  the  rays  at  the  edge  will  peel  off  and  deflect  outwards.  Thus 
the  way  that  the  beam  degenerates  is  that  the  beam  edge  gradually  erodes. 

In  this  and  the  previous  report*  we  have  studied  the  rate  at  which 

the  beam  edge  eats  its  way  in  to  the  interior.  To  do  so  we  have  considered 

the  special  case  of  an  initial  flux  distribution  having  the  shape  indicated 

in  Fig-  1.  Because  we  are  basically  interested  only  in  the  behavior  of  the 

edge  of  the  beam,  we  have  chosen  a  planar  geometry  (for  which  the  equations 
have  a  somewhat  simpler  form)  rather  than  the  more  realistic  cylindrically 
symmetrical  geometry.  The  initial  power  profile  is  taken  to  be  uniform  along 
the  negative  axis,  and  to  fall  off  exponentially  along  the  positive  axis. 

We  have  obtained  a  solution  to  the  equations  which  shows  that  the 
interior  rays  (the  rays  originating  along  the  negative  axis)  go  undeflected 
for  a  certain  distance,  then  begin  to  deflect  outwards  (towards  the  positive 
axis).  This  behavior  is  indicated  in  Fig.  2  The  solution  is  valid  only  in 
certain  regions  of  space,  but  for  any  ray  it  is  correct  up  to  and  including 
the  point  at  which  the  ray  starts  to  deflect.  Therefore,  it  provides  a  first 
order  description  of  the  desired  rate  at  which  the  beam  edge  erodes.  We  can¬ 
not  obtain  from  our  solution  the  trajectories  of  the  rays  after  they  begin 
to  peel  off. 

The  distance  from  the  laser  at  which  any  ray  begins  to  deflect  depends, 
among  other  things,  upon  its  original  distance  in  from  the  edge,  upon  the  beam 
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intensity,  and  upon  the  assured  distance-dependence  of  the  period  for  which 
any  air  mass  is  illuminated.  The  illumination  period  in  turn  depends  upon 
the  angular  velocity  at  which  the  laser  beam  is  rotated. 

For  the  model  we  lave  used,  the  degradation  of  the  beam  is  due  solely 
to  the  erosion  of  the  outer  edge.  The  interior  region  is  entirely  undegraded; 
that  is,  the  flux  in  the  interior  is  undiminished.  Without  the  introduction 
of  additional  conditions  into  the  problem  the  choice  of  a  criterion  for  when 
the  beam  has  lost  its  destructive  capabilities  is  left  somewhat  arbitrary. 

For  the  sake  of  providing  a  definite  example  we  have  required  that  the  edge 
may  be  permitted  to  eat  its  way  in  by  no  more  than  half  the  original  beam 
radius.  Using  this  condition  we  obtain  a  simple  expression  for  the  maximum 
flux  which  the  self-defocusing  effect  allows  us  to  transmit  through  the  at¬ 
mosphere.  The  expression  depends  upon  such  parameters  as  the  beam  intensity, 
the  beam  diameter,  the  target  distance,  the  rotational  rate  of  the  beam,  and 
the  effective  absorption  coefficient  for  heating  of  air  by  light  at  the  laser 
frequency 

In  summary  it  should  be  emphasized  that  our  solution  contains  a  number 
of  approximations  and  limitations,  which  have  been  mentioned  above.  Removing 
the  approximations  and  limitations  would  probably  require  an  extensive  com¬ 
puter  analysis.  The  major  value  of  our  result  relative  to  such  an  analysis 
is  that,  being  analytical,  it  provides  a  single  expression  from  which  one 
can  easily  estimate  the  maximum  transmi ttable  power  under  a  variety  of  cir¬ 
cumstances.  We  shall  ]ive  some  illustrative  examples  in  a  later  section. 
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2.  DERIVATION  OF  THE  EQUATIONS 

In  this  section  we  derive  the  equations  which  describe  the  propagation 
of  the  light  rays  in  a  medium  whose  optical  properties  are  themselves  deter¬ 
mined  by  heating  due  to  the  presence  of  the  light.  Most  of  this  derivation 
constitutes  only  a  minor  modification  of  the  work  in  reference  (1).  The 
principal  difference,  as  discussed  in  the  previous  section,  lies  in  the 
introduction  of  a  position-dependent  illumination  period  for  the  purpose  of 
simulating  the  beam  rotation. 

We  assume  a  planar  geometry  with  the  z-axis  normal  to  the  laser  face, 
the  x-axis  parallel  to  the  laser  face,  and  all  quantities  being  independent 
of  the  third  coordinate.  The  laser  face  is  in  the  z  =  0  plane,  so  that  z 
is  the  coordinate  denoting  distance  from  the  laser  along  the  propagation 
direction,  and  x  is  the  transverse  coordinate.  The  ray  trajectories  in  a 
medium  of  refractive  index  n  are  described  by  the  eikonal  equation.  One  of 
the  standard  forms  in  which  the  latter  may  be  written  is 


-here  a  is  a  unit  vector  tangent  to  the  ray  path  at  any  point,  o  is  the  path 
length  along  a  given  ray,  and  s  is  the  value  of  the  transverse  coordinate  of 
any  ray  at  z  =  0,  that  is,  x(s,z  =  0)  =  s. 

It  will  be  convenient  to  use,  instead  of  the  (x,z)  coordinate  system, 
one  defined  by  the  planes  of  constant  z  and  the  surfaces  generated  by  the  ray 
paths  (surfaces  of  constant  s),  which  are  the  solutions  to  the  problem.  Also, 
we  shall  confine  ourselves  to  situations  in  which  the  rays  deviate  from  their 
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Initial  directions  only  by  small  angles.  Therefore,  by  taking  the  projection 
of  the  vector  equation  (2.1)  along  the  transverse  axis,  transforming  to  the 
(s,z)  coordinate  system,  and  making  the  small-angle  approximation,  we  obtain 
the  ray  equation 


ax  a2x 

«  7? 


jjv; 

3s  na 

a 


(2.2) 


where  n  (s,z)  is  the  ambient  value  of  the  refractive  index. 

a 

The  conservation  of  energy  flux  F  is  stated  by  (assuming  negligible 
absorption  losses  in  the  atmosphere,  and  using  the  small-angle  approximation) 


F(s,z,t)dx  =  F(s,0,t)ds  . 


(2.3) 


This  follows  from  the  fact  that  the  surfaces  of  constant  s  are  everywhere 
tangent  to  the  ray  paths,  so  that  no  flux  escapes  through  them.  The  line 
element  dx  is  the  separation  at  height  z  between  two  rays  which  are  initially 
separated  by  ds,  so  we  may  use  the  relation  dx  =  (ax/3s)ds  to  obtain  from 
(2.3)  the  flux  conservation  equation, 

*  <2-4> 

Our  neglect  of  absorption  losses  is  quite  justifiable  for  ruby  and  neodymium 
lasers,  for  which  the  energy  loss  in  passing  vertically  through  the  entire 
atmosphere  is  no  more  than  about  one  percent.  For  a  COj  laser  about  0.3  of 
the  incident  energy  is  absorbed  in  the  atmosphere.  This  loss  is  not  insigni- 
cant,  but  for  our  purposes  it  is  not  great  enough  to  warrant  the  additional 
complications  of  accounting  for  it  in  the  conservation  equation. 


Next  we  need  the  relation  between  the  flux  and  the  refractive  index. 
The  rate  of  change  of  refractive  index  Is  proportional  to  the  rate  of  change 
of  the  density,  which  is  proportional  to  the  rate  at  which  heat  Is  deposited 
in  the  atmosphere;  the  latter  In  turn  is  proportional  to  the  flux.  The 
refractive  index  for  a  dilute  gas  is  given  by 


n-1  =  2iTpa 


(2.5) 


where  p  is  the  molecular  number  density  and  a  is  the  molecular  polarizability. 
Consequently,  we  have 


dn  0  do  n-1  dp 
ot  ’  2™Ht  =  ~T  it 


(2.6) 


The  density  changes  according  to  the  equation  of  state  of  an  ideal  gas, 

P V  =  NKT  ,  (2.7) 


where  p  =  N/V,  and  where  K  denotes  Boltzmann's  constant.  We  shall  assume 
that  the  heating  takes  place  at  constant  pressure,  so  that  from  (2.7)  we 
obtain 


1  dp  1  dV  NK  dT 

p  Ht  “  T  cTt  '  “  W  31  ‘ 


(2.8) 


The  assumption  of  isobaric  heating 
interest  ere  large  compared  with  the  time 
heated  region.*  The  velocity  of  sound  in 


is  valid  only  if  the  times  of 
required  for  sound  to  traverse  th* 

4 

air  is  about  3x10  cm/sec,  so  it 
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would  take  several  milliseconds  to  cross  a  beam  a  meter  in  diameter.  This 
means  that  If  we  are  considering  a  stationary  target  the  pulse  length  must  be 
greater  than  several  milliseconds,  or  less  for  a  smaller  beam.  For  a  moving 
target  we  have  the  additional  constraint  that  the  beam  must  rotate  sufficiently 
slowly  that  It  Is  effectively  stationary  during  a  sound  traversal  time.  It 
may  easily  be  seen  that  this  constraint  is  equivalent  to  the  requirement  that 
the  target  velocity  be  somewhat  less  than  the  velocity  of  sound  In  air. 
Therefore  our  results  will  be  applicable  only  to  the  illumination  of  subsonic 
targets.  The  direction  in  which  our  results  err  for  supersonic  targets  is 
easy  to  see.  When  a  region  of  atmosphere  is  heated,  a  finite  amount  of  time 
is  required  for  the  density  change  to  develop.  This  is  the  time  needed  for 
sound  to  cross  the  heated  region.  Over  a  shorter  period  the  density  change 
and  therefore  the  refractive  index  change, ic  less  than  what  we  estimate 
under  the  assumption  of  heating  at  constant  pressure.  Consequently  the  de- 
focusing  effect  will  be  smaller  than  estimated.  As  a  result,  the  amount  of 
flux  that  can  be  transmitted  in  the  case  of  a  supersonic  target  will  be 
greater  than  the  upper  limit  which  we  shall  predict. 

To  continue,  the  temperature  increase  dT  for  an  amount  of  heat  dQ 
deposited  at  constant  pressure  in  a  volume  V  of  air  is 

dQ  =  CpdT  ,  (2.9) 

where  Cp  Is  the  specific  heat  at  constant  pressure.  Now  if  we  denote  by  k  the 
effective  absorption  coefficient  for  heating  of  air  by  the  laser  light,  then 
the  amount  of  heat  per  unit  volume  per  unit  time  deposited  due  to  beam  absorp¬ 
tion  is  given  by  kF,  Therefore,  from  (2.9)  we  obtain 
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(2.10) 


CP  dT 

=  r  at 

Now,  combining  (2.10),  (2.8)  and  (2.6),  and  using  the  well  known  ideal  gas 
relation, 


NK  = 


(2.11) 


Cy  being  the  specific  heat  at  constant  volume,  we  obtain  the  desired  relation 
between  refractive  index  and  flux, 


dn 

Ht 


(n-1) (y-l)kF 
yP  * 


(2.12) 


where  y  =  Cp/Cy.  For  air  the  ratio  of  specific  heats  is  y  =  1.4.  Eq . ( 2 . 1 2 ) 
may  be  Integrated  to  give 

F  [(s(x,z,t)  ,2,t]dt  ,  (2.13) 

x,z 


n  ( s ,  z .  t )  -  n 


where  n,  is  the  ambient  value  of  the  refractive  index  in  the  absence  of  the 
beam.  The  factor  n-1  on  the  right  hand  side  of  (2.12)  has  been  assigned  the 
fixed  value  n-1  and  removed  from  the  integral,  since,  as  may  be  easily 

a 

verified,  the  fractional  change  in  n  is  negligible  compared  to  unity.  The 
subscript  x,z  on  the  integral  sign  emphasizes  that  the  relation  (2.12) 
applies  to  a  fixed  air  mass,  and  therefore  to  a  fixed  point  in  space  (x,z). 
Since,  because  of  continued  heating  of  the  air,  the  rey  trajectories  change 
with  time,  the  origin  coordinate  s(x,z,t)  of  the  ray  passing  through  the 
point  (x,z)  also  changes  with  time.  Thus,  in  the  integral  f  depends  on  t  not 
only  explicitly,  but  also  implicitly  through  the  implied  variation  of  s  with 
time. 
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Now  we  make  our  most  serious  approximation.  We  assume  that  in  the 
Integral,  F  may  be  regarded  as  constant  with  respect  to  both  t  and  s  during 
the  period  for  which  the  beam  passes  through  the  given  air  mass.  We  replace 
the  integral  by 

F[s(x,z,t),z,t]dt  F(s,z,t)r(z,t)  ,  (2.14) 

'X,Z 

where  r(z,t)  is  the  time  interval  over  which  an  air  mass  at  distance  z  from 
the  laser  has  been  illuminated  at  time  t.  How  well  we  simulate  the  rotation 
of  the  laser  beam  now  depends  on  how  judiciously  we  choose  r(z,t).  In  the 
special  case  of  a  stationary  target  we  have  i(z,t)  *  t. 

Strictly  speaking,  the  beam  rotation  should  be  described  by  time 
dependent  boundary  conditions  on  the  rays.  What  we  shall  do  here,  as  in  th«* 
stationary-target  case.  Is  to  let  the  z  *  0  plane  represent  the  laser  face 
and  require  th:t  the  rays  emerge  normally  from  the  laser  face.  That  is,  the 
angle  made  by  the  rays  relative  to  the  ncrr?.1 ,  -  =  ^x/'-z,  must  be  zero  at 
z  =  0  To  represent  a  rotating  beam  properly  we  must  do  something  much  more 
complicated  mathematically.  That  is,  we  must  require  that  the  rays  emerge 
normally  from  a  plane  which  is  rotating  relative  to  a  coordinate  system  fixed 
in  the  atmosphere.  Our  use  of  the  approximation  (2.14)  eliminates  this  com- 
pl ication 

Now,  according  to  (Z.5)  n  -1  is  proportional  to  .  ,  the  atmospheric 

u 

density,  and  c  decreases  exponentially  with  altitude.  Mowever,  according  to 
the  ideal  gas  law  (2.7)  the  pressure  P  is  'Iso  proportional  to  *.  Consequently, 
in  (2.13)  the  ratio  (n  -1)/P  is  independent  of  density,  and  therefore  of 

a 
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altitude,  so  that  the  sea  level  values  of  nfl  and  P  may  be  used.  In  general 
the  altitude  dependence  of  k,  the  effective  absorption  coefficient  for  heating, 
Is  considerably  more  complicated  and  uncertain.  For  example,  for  ruby  lasers 
essentially  all  of  the  absorption  is  by  water  vapor,  and  for  CO2  lasers  about 
half  Is  by  water  vapor.  The  distribution  of  atmospheric  water  vapor,  however, 
varies  drastically  from  one  geographical  location  to  another,  as  well  as  from 
one  time  to  another.  Therefore,  for  the  moment  we  shall  simply  regard  k  as 
an  unspecified  function  of  altitude,  that  Is, 

k  =  k(zcos$)  ,  (2.15) 

where  #  Is  the  zenith  angle  at  which  the  beam  is  propagating  through  the 
atmosphere  at  any  instant  of  time. 

We  now  have  for  the  relation  between  the  refractive  Index  and  the  flux, 

(>-l)(n  -l)k(7cos$)F(s,z,tMz,t) 

n(s.z.t)  -  na(z)  = - - - ^ -  .  (2  16) 


Now,  the  system  of  equations  which  we  wish  to  solve  is  composed  of  the 
ray  equation  (2.2),  the  flux  conservation  equation  (2.4),  and  the  relation 
between  flux  and  refractive  index  (2.16).  First  we  eliminate  the  flux  between 
(2.4)  and  (2.16),  and  obtain  a  relationship  between  the  value  of  the  refractive 
index  at  any  point  (s,z)  and  its  value  along  the  same  ray  path  but  at  the 
lase’-  face, 


n(s ,Z ,t )  -  n ^(z) 


n($,C)-n  (0) 

_ _ d _ 

•x/  s 


TU.t)  . 


(2. i 7 ) 


64 


where  wp  have  defined 


T(z,t) 


k(zco.,i)  i(z,t) 

~T~0  TltSVtT  * 


(2.18) 


and  k  k(0). 
o 

Next,  it  is  convenient  to  define  the  dimensionless,  positive  quantity, 


f(s,t)  = 


n  (0)-n(s,0,t) 

“srnn 


(2.19) 


so  that  (2.17)  may  be  rewritten  in  the  form  ^recall  that  n  is  so  close  to  unity 
that  unless  index  differences  are  involved  we  may  write  ng(z)  =  na(0)] 


na(z)-n(s,z,t) 


f(s  ,t)T(z,t) 
3x/3S 


(2.20) 


Now,  physically,  f(s,t)  may  be  regarded  as  a  reduced  flux.  To  see  this  we 
use  (2  1 G )  and  the  definition  (2.19)  of  f(s,t)  and  obtain 


f(s,t) 


(-,-l)(na-l)koT(0,t) 

- _p 


F(s,0,t)  . 


(2.21) 


rhus,  f(s,t)  is  seen  to  have  the  same  s-dependence  as  F(s,0,t).  That  is,  it 
is  proportional  to  the  power  profile  at  the  laser  face. 

Finally,  we  may  eliminate  the  refractive  index  by  inserting  (2.20)  into 
the  ray  equation  (2.2)  to  obtain 


3X  3^X  1 

a 

[f(s)  ' 

as  dz?.  T(zT 

'  85 

ax/as 

(2.22) 
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We  have  simplified  the  notation  by  suppressing  the  parametric  t-dependence. 
To  obtain  the  ray  trajectories  we  must  now  solve  (2.22)  with  the  appropriate 
boundary  conditions. 


x(s,G)  =  s  ,  (2.23a) 

e(s,0)  =  !f(s,0)  =  °*  (2.23b) 
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3.  SOLUTION 


Performing  the  derivatives  in  (2.22)  and  setting 


X  =  S  +  Xj  , 


(3.1) 


(so  that  x^  is  the  transverse  displacement  of  the  ray  originating  at  x  =  s) 
we  may  rewrite  (2.22)  in  the  forr 


1  ♦  T 


2?\  1  f(  ,  :'2>(1  df /,  .  3xli  n 

t-cft  ' f(s)  ~7“  a?!1  *  it)  "  0 


(3.2) 


To  solv?  (3.2)  we  shall  make  use  of  the  following  transformations  of  variables, 


/f(s) 


§  ■ 


(3.3) 


(3.4) 


Actually,  it  will  be  useful  to  transform  only  selected  terms  in  (3.2),  so  as 


to  recast  it  in  the  form 


sxi|3|  Ai 


IV  1  ii  Il'.l 
3^  2T7  32  32  I 


L  3  3xl]  df 

^r+  y  +  7  arj h? 


(3.5) 


where 


Xj(p»0  =  x^(s,z) 


(3.6) 


3X,  .  y  9X-. 

We  now  assume  that  in  (3.5)  the  — —  terms  and  the  -=-*  —  t— -  term 

u  S  2J6  v  Z  0  Z 

may  be  dropped.  Both  of  these  terms  vanish  at  the  boundary  plane  z  *  0, 
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ax, 

because  of  the  boundary  conditions  Xj(s,0)  =  0  and  e ( s ,0)  =  (s,0)  =  0. 

We  shall  see  that  in  fact  there  exists  a  useful  solution  for  which  they  vanish 
over  an  entire  region  of  space;  in  that  region  of  space  the  solution  will  be 
exact. 

The  approximate  equation  which  we  obtain  from  (3.5)  is  now 


A  _  iV  df 


=  0 


95 


9p 


(3.7) 


Equation  (3.7)  is  linear  with  constant  coefficients,  and  can  be  solved  in 
closed  form.  We  shall  net  go  through  the  details  of  solving  it  here;  however, 
for  completeness  the  procedure  has  been  included  in  an  appendix.  The  solution, 
with  the  required  boundary  conditions,  is 


x^ (s ,z)  =  2s  -  g(s,s)  -  g(s,-5)  , 


(3.8) 


where  g(s,0  is  given  by 


g(s,c) 
ds ' 
/f(s') 


(3.9) 


and  5  ir  turr.  i:  a  "unction  of  z  according  to 


5(z)  =  ^  A< z' )  dz1  . 


(3.10) 


The  last  equation  follows  from  ..he  transformation  (3.4) 
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Now,  the  terms  we  have  dropped  do  not  vanish  in  general.  It  Is  only 
for  a  special  choice  of  the  initial  power  profile  f (s )  that  they  vanish  In  a 
region  of  interest.  That  special  case  was  mentioned  in  the  introduction  and 
is  illustrated  in  Fig.  1.  The  initial  flux  distribution  is  uniform  along  the 
negative  axis,  s  <  0,  which  represents  the  interior  of  the  more  properly  cy¬ 
lindrical  beam.  It  falls  off  exponentially  along  the  positive  asix;  the  fall- 
off  region  represents  the  beam  edge.  Specifically,  we  have 


f(s) 


s  «  0 


e-2s/x 


s  *  0 


(3.11) 


for  times  greater  than  t  =  0.  For  t  <  0,  we  assume  the  flux  to  be  zero;  that 

p 

is,  the  beam  is  turned  on  at  t  =  0.  The  constant  a  ,  according  to  (2.21),  is 
given  by 


2  (Y-l)(na-l)koFoT(0,t) 
'  =  yP 


(3.12) 


where  FQ  is  the  (uniform)  initial  flux  magnitude  along  the  negative  transverse 
axis. 

It  is  easy  to  see  physically  why  we  expect  that  for  this  flux  distri¬ 
bution  the  terms  we  have  dropped  from  (3.5)  will  vanish  over  a  finite  region 
of  space.  Along  the  negative  axis  the  flux  is  initially  uniform,  so  the  air 
will  be  heated  uniformly.  Consequently,  the  refractive  index  will  change,  # 

but  its  change  will  be  uniform  with  respect  to  the  transverse  direction,  so 
there  will  be  no  transverse  gradient.  As  a  result,  the  light  rays  will 
initially  not  be  bent,  so  the  displacement  x^  will  be  zero;  furthermore,  the 
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uniformity  will  propagate  in  the  z-direction.  In  the  region  where  has  the 
constant  value  zero,  we  have  ax^/as  =  ax^/az  =  0,  so  that  the  terms  we  have 
ignored  actually  vanish  identically. 

The  situation  is  different  along  the  positive  transverse  axis.  The 
flux,  and  therefore  the  heating  rate,  is  greatest  at  the  beam  edge  (s  =  0) 
and  falls  off  for  larger  positive  s.  Therefore,  the  refractive  index  increases 
in  the  positive  transverse  direction.  According  to  the  eikonal  equation  (2.1), 
the  rays  deflect  in  the  direction  of  increasing  index.  Therefore,  they  im¬ 
mediately  begin  to  peel  off  in  the  positive  transverse  direction,  which  is 
equivalent  to  peeling  off  outward  from  a  cylindrically  symmetrical  beam.  As 
these  outer  rays  peel  off,  the  beam  edge  is  displaced  inwards  (in  the  negative 
transverse  direction)  with  increasing  z,  and  the  rays  at  the  new  edge  peel 
off.  Thus  the  rays  in  the  region  of  initially  uniform  flux  distribution 
propagate  undeflected  only  for  a  limited  distance,  then  they  deflect  outwards 
(in  the  positive  transverse  direction). 

To  see  just  how  far  they  propagate  undeflected  we  insert  the  explicit 
expression  (3.11)  for  the  initial  flux  distribution  into  the  general  solution 
(3.8,  3.9,  3.10)  and  obtain  the  explicit  solution, 


f  0 


s  *  -as 


x^s.z)  = 


s+as-Aloge(l+s/A+as/x)  , 
2s*A+as->es^x-Aloge(as/x+es^x) 


l  2s-Aloge(e2s/A-aV/A2]  , 


-as  <  s  <  0 
0  <  s  «  Alogg( 1+as/x ) 
*1oge(l+as/A)  <-  s  .  (3.13) 
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Here  f(z)  is  understood  to  be  a  known  function  of  z,  determined  from  Eq.  (3.10). 
To  obtain  an  explicit  expression  for  ^(z)  we  must  first  specify  T(z,t);  doing 
so  is  equivalent  to  specifying  the  mean  illumination  time  of  an  air  mass  at 
a  distance  z  at  time  t.  In  order  to  get  a  feel  for  the  properties  of  the 
solution  (3.13)  it  is  probably  best  to  think  of  it  in  terms  of  the  simple 
special  case  of  a  stationary  target,  with  both  the  target  and  the  laser  at 
zero  altitude  In  this  case  we  find  from  (2  18)  that  T(z,t)  =  1,  and  as  a 

result  we  have  %{z)  -  z 

In  the  'egion  s  <•  -a l(z),  then,  we  have  3Xj/3z  =•  sx^/ss  =  0,  so  that 
in  this  region,  as  we l ’  as  everywhere  along  the  boundary  plane  z  =  0,  (3.13) 
is  an  exact  solution  of  (35)  Elsewhere,  as  we  have  said,  the  rays  deflect, 
so  that  the  terms  we  have  dropped  in  going  from  (3.5)  to  (3.7)  do  not  vanish. 

In  fact,  these  terms  quickly  become  sufficiently  large  that  the  solution  (3.8) 
is  not  a  good  approximation  except  in  and  very  close  to  the  region  s  -a«;(z). 
For  s  -aeU).  but  at(z)  -  jsj  ^  l,  it  may  be  seen  that  the  neglected  terms 
a^e  st,n  quite  smai',  so  that  the  solution  is  still  quite  accurate.  From 
(3  13)  we  see  that  under  these  conditions  the  ray  trajectories  are 

*j(s,z)  r  75--  (ac+s)^  -  (ae - 1 s |  (3.14) 

Therefore,  we  find  that  the  outer  edge  of  the  beam  (the  region  of 
diverging  rays)  eats  into  the  interior  region  (the  region  of  initially  uniform 
flux  distribution)  along  the  curve  s  -  -ac(z)  That  is,  for  any  negative  s, 
the  corresponding  ray  propagates  undeflected  until  it  reaches  a  distance  z 
from  the  laser  face  given  by  %(z)  -  -S/a;  then  it  begins  to  diverge  para- 
bolicaily  m  e.(z)  according  to  13  14)  This  behavior  is  illustrated  in  Fig. 
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It  should  be  noted  that  the  distance  z  at  which  an  interior  ray  first  begins 
to  deflect  does  not  depend  on  the  steepness  of  the  beam's  original  edge  shape. 
However,  once  the  ray  begins  to  deflect,  it  does  so  more  rapidly  for  a  more 
steeply  descending  beam  edge;  more  precisely,  we  see  from  (3.14)  that  the  ray 
displacement  x^  is  proportional  to  1/x,  where  \  is  the  scale  length  for  the 
exponential  fall -off  of  the  beam  edge. 

The  region  in  which  (3.14)  is  valid  is  quite  small;  from  (3.14)  we  only 

learn  the  initial  behavior  of  the  rays  when  they  begin  to  deflect.  Since  we 

do  not  know  the  details  of  how  much  the  beam  spreads  once  it  does  so,  we  must 

make  a  somewhat  arbitrary  choice  of  a  criterion  for  when  the  beam  is  to  be 

considered  too  degraded  to  be  effective.  Now,  we  have  seen  from  the  solution 

(3.13)  that  the  beam  interior  is  entirely  undegraded,  in  the  sense  that  the 

power  density  is  undiminished;  however,  the  beam  edge  moves  steadily  inward, 

so  that  the  undegraded  interior  decreases  in  size.  We  assume  that  when  the 

region  of  undiminished  power  density  gets  too  small  the  beam  is  no  longer 

useful.  Therefore,  if  we  denote  by  s  the  maximum  distance  which  we  can 

o 

tolerate  having  the  beam's  edge  eat  into  its  interior,  then  from 

S(z)  =  sQ/a  (3.15) 

and  the  definition  (3.12)  of  a,  we  find 

F°  ’  (v-l)(na-l)ko,(0.t){Z(z)  '  U  '6) 

For  a  specified  sQ  and  a  target  at  distance  z,  FQ  is  then  the  limit  due  to 
self-defocusing  on  the  flux  which  can  be  transmitted  through  the  atmosphere 
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It  is  desirable  to  simplify  (3.16)  somewhat  by  going  back  to  the 
definitions  (3.10)  of  c(z)  and  (2.18)  of  T(z,t)  and  noting  that 


C(z)  = 


1 


AQT(0,t) 


■z  ^ 

[k(z 'cos4>)t (z*  ,t)]  dz '  . 


(3.17) 


The  term  kQT(0,t)  in  (3.16)  may  be  canceled  and  we  obtain  finally 


ypso 

(y"l)(n  -l)n^(z) 


(3.18a) 


where 


n(z)  = 


[k(z ' cos4> )t (z *  ,t )]  dz* 


(3.18b) 


Before  we  can  obtain  numbers  from  (3.18)  we  must  still  choose  a  function 
x(z,t)  which  describes  the  beam  rotation,  and  we  must  pick  a  value  of  sQ.  In 
the  next  section  we  shall  give  several  illustrative  examples  of  the  use  of 
(3.18).  In  all  of  these  examples  we  shall  assume  that  the  beam  has  lost  its 
destructive  ability  when  the  beam  edge  has  eaten  its  way  in  by  one  quarter  of 
its  original  diameter,  sQ  =  D/4. 


f 
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4. 


ILLUSTRATIVE  EXAMPLES 


In  this  section  we  shall  treat  several  special  cases  so  as  to  illustrate 

the  use  of  (3.18)  In  calculating  FQ,  the  flux  limit  due  to  self-defocusing. 

* 

We  shall  assume  we  are  using  a  COg  cw  laser.  We  mentioned  in  Section  2  that 
the  altitude  dependence  of  the  absorption  coefficient  for  C02  Is  rather  com¬ 
plicated  and  uncertain.  We  shall  simply  assume  that  for  our  purposes  it  may 
be  adequately  represented  by  the  altitude  dependence  of  the  atmospheric 
density, 

k(zcos$)  =  kQ  e-(z/L}cos*  ^  (4.1) 

where  L  is  the  atmospheric  scale  height,  L  =  8  km;  the  value  of  kQ  is 
kQ  =  3xl0”7cnf ‘.  The  values  of  the  other  constants  which  appear  in  (3.18) 
are  as  follows:  The  ratio  of  specific  heats  in  air  is  \  -  1.4,  the  atmospheric 
pressure  is  P  *  10  dynes /cm  ,  and  the  ambient  refractive  index  is 
n  -1  •-*  3xl0~*.  Also,  In  our  examples  we  shall  take  the  Initial  beam  diameter 

O 

to  be  0  *  lm 

Case  1.  Stationary  Target 

For  a  stationary  target  eacn  air  mass  is  illuminated  for  the  entire  laser 
pulse  length  t,  so  that  we  have 


*  It  is  worth  recalling  here  that,  as  was  pointed  out  in  Section  2,  there  *s 
essentially  no  self-defocusing  problem  for  giant  pulse  lasers;  the  density 
change  requires  a  finite  amount  of  time,  which  is  of  the  order  of  milliseconds 
for  a  beam  of  1  meter  diameter.  The  current  giant  pulse  laser  is,  however,  not 
the  answer  to  the  self-defocusing  problem,  since  the  amount  of  energy  »hat  can 
be  put  out  in  one  pulse  is  In  general  far  too  small  for  most  applications.  If 
we  think  in  terms  of  a  rapid  succession  of  giant  pulses,  then  we  are  in  effect 
back  to  the  cw  laser. 


74 


T (Z , t)  =  t 


(4.2) 


Using  (4.2)  and  (4.1)  to  obtain  n(z)»  and  taking  sQ  =  D/4,  where  D  is  the 
beam's  initial  diameter,  we  immediately  find  from  (3.18), 


F  t  s  _ yPO^cos2^ _ 

0  64(y-l)(n  -l)knL2(l-exp(-zcos^/2L)? 


(4.3) 


Note  that  the  limitation  here  is  not  on  the  flux  F  ,  but  rather  on  the  total 
energy  per  unit  area  Ft  which  can  be  transmitted. 

If  we  consider  the  case  of  a  target  in  the  zenith  direction  U  *  0) 
at  an  altitude  much  greater  than  t.he  atmospheric  scale  height  (z  >>  L),then 
(4.3)  simplifies  to 


Ft 


,PD£ 


54(v-l)(iyl)k0L‘ 


(4.4) 


Putting  in  the  numbers  specified  above  we  obtain 

Ft  8  .95  joules/cm  . 

Now  let  us  consider  the  case  in  which  both  the  target  and  the  laser  are 
on  the  ground,  so  that  coss  *  0.  In  this  case  (4.3)  reduces  to 


V 


JL 


? 

PD“ 


(4.5) 


For  a  target  at  a  distance  of  z  *  1  km,  (4.5)  yields 


Ft  *  24C  joules/on  . 
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Case  2.  Moving  Target 

For  a  moving  target  the  accuracy  of  our  results  depends  on  how  judi¬ 
ciously  we  choose  the  expression  for  t(z,t),  the  mean  illumination  time  of 
an  air  mass  at  distance  z  from  the  laser.  If  the  angular  coordinate  0  of  the 
target  Is  changing  at  the  rate  9,  then  at  large  distances  z  from  the  laser 
the  beam  is  moving  with  a  linear  velocity  ez,  so  that  the  time  taken  for  the 
beam  to  cross  the  air  mass  is 


t(z,t)  -  D/ez  . 


(4.6) 


At  snail  distances  this  expression  is  quite  invalid.  The  concept  of  an 
average  illumination  time  is  meaningful  only  if  this  time  is  less  than  the 
total  pulse  time,  that  is,  if  i(z,t)  <  t.  However,  the  above  expression  for 
t  blows  up  at  small  z,  so  the  required  inequality  is  clearly  violated  If, 
nevertheless,  we  Insert  (4.6)  into  the  expression  (3.18b)  for  n(z)  to  obtain 


r(z) 


S 

k  D  \  rZ 

-O-l  ! 


e-(z'/2L)cos* 

/? 


dz* 


(4.7) 


we  may  note  that  for  small  z'  the  integral  gees  like  *  z ' ;  consequently,  the 
contribution  to  the  integral  is  small  from  the  range  in  which  (4.6)  over¬ 
estimates  the  average  Illumination  time. 

We  shall  use  (4.7),  but  first  let  us  digress  momentarily  to  get  a 
crude,  but  somewhat  more  quantitative  estimate  of  the  error  involved.  For 
simplicity  we  take  cos#  1  0,  so  that  the  exponential  in  the  integral  becomes 
unity  and  the  integral  is  directly  evaluable.  Then  if  we  regard  (4.6)  as 
invalid  up  to  a  distance  zQ,  but  va*id  for  all  larger  distances,  the 
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fractional  error  is  immediately  seen  to  be  c 
at  which  i  becomes  as  large  as  t,  so  that  frc 
However,  from  (4.6)  we  also  find  that  0/*  =  2 
error  may  be  estimated  to  be  v'z0'z  =  *  t  ( 2 » t )/ 


N  '•>,  ;  is  the  distance 


.6)  *•  ive  z  -  D/-t. 

0 

t).  -nefore,  the  trsct 
,ie  co'  ude  then,  that  i# 


the  total  pulse  time  t  is  sufficiently  greate  ^  tn*  veraqe  Illuminati 


time  :(z,t)  at  the  target,  we  may  use  (4.7). 
Returning  to  (4.7)  we  obtain,  after  a  s 


tra'  ormatlon  on  the 


integral , 


where 


/  2'Lk  D\  _ 
n(z)  =  I— -  E("ZCO 

\  acos;/ 


E(X)  -  4=  f  e'^2 

/2r  i-X 


is  the  normal  probability  integral,  having  th«  roperty  (®)  -  1.  Now  from 
(3.18a)  we  obtain  the  desired  expression  for  1  limits  ux. 


.  P0-.-COS; 


32-( ,-l)(na-llkoLE<:(.zcc 


If  the  target  is  flying  directly  overhe* 
so  that  (4.10)  becomes 


■tave  s:  5  1  and  ••  *  v/ 


32-(.-l)(na-l)koLz  £2[  .z 
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Let  us  assume  that  the  target  is  at  an  altitude  z  =  10  km  and  has  a  velocity 
v  =  500  mph  =  2.23  x  10^cm/sec.  We  have  assumed  a  subsonic  velocity  (sound 
velocity  in  air  is  3  x  10^cm/sec),  since,  as  we  discussed  earlier,  our  model 
is  not  applicable  to  supersonic  velocities.  For  the  chosen  altitude  and 
velocity  we  find  from  (4.11) 

Fq  =  197  watts/cm2. 

The  average  illumination  time  of  an  air  mass  in  the  neighborhood  of  the  target 

_3 

is  t  =  D/v  =  4.5  x  10  sec,  so  this  result  for  FQ  is  meaningful  only  if  the 
total  laser  pulse  time  length  t  is  considerably  greater  than  several  millisec. 
For  laser  pulses  shorter  than  x  =  D/  one  would  be  more  accurate  in  regarding 
the  beam  as  stationary  and  using  (4.3). 

If  the  target  is  on  the  ground  and  moving  normal  to  the  line  of  sight 
with  velocity  v,  we  again  have  e  =  v/z,  but  this  time  cos<t>  =  0,  so  that  (4.10) 
reduces  to  (since  E(x)  -  ^2/7  x  for  x  <<  1) 


F 


o 


yPDv _ 

64(Y-l)(na-l)koz2 


(4.12) 


Taking  a  target  velocity  v  =  10  m/sec  (=  21.6  mph)  and  a  target  distance 
z  =  1  km,  we  find 

Fq  =  610  watts/crn2. 


In  this  case 
is  t  =  D/v  = 
greater  than 


the  average  illumination  time  of  an  air  mass  near  the  target 
10’1  sec,  so  that  the  laser  pulse  length  must  be  somewhat 
10'1  sec  for  our  model  to  apply. 
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APPENDIX 


The  equation  to  be  solved  is 


3  X1  3  X1  df  n 

— 2 - r  Hs  ’  0  ' 

H  ap 


(Al) 


where  Xj  (p,c)  is  the  ray  displacement, 

X^p.5)  =  x^s.z) 


(A2) 


and  the(p,c)  coordinate  system  is  related  to  the  (s,z)  coordinate  system  by 
the  transformations 


=  1 
/fU) 

(A3) 

§ 

=  /r(?) 

(A4) 

The  boundary  conditions  are  that  the 

initial  displacement  is  zero, 

Xj (s ,0) 

=  o 

(A5) 

and  the  initial  angle  of  any  ray  relative  to  the  z  axis  is  zero, 

9(s,0) 

•  ^  (s  ,0)  =  0  . 

(A6) 
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Now. we  define  the  functions  S(p),  P(s)  and  c(z)  by 


fS(p)  ds ' 
o  /?(F) 

(A7) 

p  ■  J 

P(s)  =  1 

fs  ds1 

0  v'fTs’l  • 

a  _ _ _ 

(A8) 

C(z)  = 

/T(z‘ )  dz ’ 

0 

(A9) 

According  to  (A3)  the  s-p  transformation  and  its  inverse  are  then  given  by 


s  =  S(p)  ,  p  =  P(s)  .  (A10) 

It  will  be  seen  shortly  that  defining  the  functions  S(p)  and  P(s)  is  useful 
for  avoiding  notational  confusion. 

Next  we  transform  the  last  term  in  (Al)  as  follows,  using  (A3)  and 

(A10), 


df 

a? 


f  * «S(P)1  ■  ^ 


d/F 


,  d  ds  ,  d2S 
2  Hp  Hp  =  2  “ 


dp 


ifi 


(All) 


so  that  (Al)  becomes 


?  d2S 

7T“  7T"  +  2  7T 

dE,  9p  dp 


=  0 


(A12) 
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The  general  solution  of  (A13)  is  well  known  to  be 


XjCp.c)  -  2S(p)  «  A(p+c)  +  B(p-c)  ,  (A14] 

where  A  and  B  are  arbitrary  functions. 

Now  we  evaluate  A  and  B  from  the  boundary  conditions  (A5)  and  (A6). 
Since  at  z  =  0  we  have  5  =  0,  it  follows  from  (A5)  that 

Xj (p ,0)  =  2S(p)  -  A(p)  +  B(p)  =  0  ,  (A15) 

so  that 

B(p)  =  -2S(p)  -  A(p)  .  (A16) 


We  may  now  write  (A14)  as 


X^P.C)  =  2S(p)  +  A(p+c)  -  2S(p-^)  -  A ( p- ^ )  .  (A17) 

To  apply  the  second  boundary  condition  (A6)  we  use  the  transformation 

(A4), 

3X  9X 

*r  *  w  *  A'(p+0  +  A'(p-e)  +2S'(p-0  .  (A18] 

Then  from  (A6)  we  obtain 

2A'(p)  +  2S'(p)  =  0  ,  (A19! 
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or 


A(p)  =  -S(p)  +  c  ,  (A20) 

where  c  Is  an  arbitrary  constant.  Inserting  this  result  into  (A15)  we  now 
have 

Xj(p,0  =  2S(p)  -  S(p+5)-S(p-0  ,  (A21) 


The  final  task  is  to  transform  back  to  the  (s,z)  coordinate  system. 

From  the  transformations  ( A7 ,8 ,9 )  we  have 

Xj  (s  ,z)  =  2s-StP(s)n(z)]-S[P(s)-5(z)]  .  (A22) 

In  addition  we  have 


Subtracting 


P(s)H(z) 

P(s)H(z) 

(A24)  from  (A23) 

5(z) 


ds* 


/W) 


+  S(z)  , 


fSlP(s)+c(z)l  ds 
0 


f(s') 


yields 


■f 


stP(s)n(z)l  ds, 


(A23) 


( A24 ) 


(A25) 


which  directly  determines  the  function  SlP(s)+c(z)l.  Finally,  we  simplify 
the  notation  by  defining  the  function  gls,c(z))  as  follows, 

g(s.c)  =  S[P(s)H(z)l  .  (A26) 
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The  solution  may  now  be  written  In  the  (s,z)  coordinate  system  as 


x^s.z)  =  2s-g(s,5)-g(s-5)  . 


(A27) 


where  g(s,  )  Is  determined  from 


g(s  ,c;  ds. 

S 


(A28) 


and  c(z)  Is  defined  as  a  function  of  z  by 


,z  _ 

?(z)  =  A(z')  dz* 


(A29) 
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Fig.  1.  Uniform  distribution  of  flux  FQ(s)  emanating  from  a  semi¬ 
infinite  laser  face  (s  <  0),  with  an  exponential  fall-off 
at  the  edge.  The  rays  originating  in  the  uniform  region 
along  the  negative  s-axls  are  initially  undeflected. 

Those  originating  along  the  positive  s-axls  are  iamediately 
deflected  to  the  right. 


Fig.  2.  Trajectories  of  rays  emanating  from  the  region  of 
uniform  flux  distribution  (s  <  0)  depicted  In 
Fig.  1.  The  rays  go  out  vertically  until  they 
Intersect  the  dashed  line,  then  they  begin 
deflecting  to  the  right  parabollcally  with 
distance  z. 
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in  the  case  of  high  energy  laser  beam  propagation  in  the  atmosphere.  Two  subjects  are 
investigated,  namely,  heating  mechanisms,  and  the  dynamics  of  the  propagating  laser 
beam.  The  heating  mechanism  study  is  concerned  specifically  with  photo-absorption  in 
the  far  wings  of  pressure  broadened  CO?  vibration-rotation  bands,  while  the  beam  dy¬ 
namics  analysis  treats  the  effects  of  target  motion  on  the  degree  of  self-defocusing 
of  the  laser  beam.  A  theory  of  far-wing  pressure  broadening  is  developed  in  which  the 
basic  broadening  mechanism  is  assumed  to  arise  from  perturbations  of  the  absorber's 
rotational  motion  via  a  short  range  repulsive  interaction  with  a  colliding  molecule. 
The  theory  predicts  a  far-wing  spectral  behavior  described  by  a  product  of  a  lorentz 
line  shape  and  an  exponentially  decreasing  factor,  in  general  agreement  with  recent 
measurements.  In  the  beam  dynamics  work  a  convenient  closed  form  expression  has  been 
obtained  for  the  maximum  flux  which  the  thermal  self-defocusing  effect  permits  one  to 
transmit  through  the  atmosphere.  This  analysis  takes  into  account  the  effect  of  tar¬ 
get  motion.  The  expression  for  the  maximum  flux  depends  upon  such  parameters  as  the 
initial  beam  intensity,  the  initial  beam  diameter,  the  target  distance,  the  rotational 
rate  of  the  beam,  and  the  effective  absorption  coefficient  for  heating  of  air  by  light 
at  the  laser  frequency. 
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